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O Abstract 

<^ ' We systematically analyze the supersymmetric contributions to the mixing CP 

asymmetries and branching ratios of B — > <j)Kg and B — ► rj ' K$ processes. We 
consider both gluino and chargino exchanges in a model independent way by using 
P-i. the mass insertion approximation method. While we adopt the QCD factoriza- 

tion approach for evaluating the corresponding hadronic matrix elements, a critical 
comparison with predictions in naive factorization one is also provided. We find 
that pure chargino contributions cannot accommodate the current experimental re- 
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suits on CP asymmetries, mainly due to b — > constraints. We show that charged 
Higgs contributions can relax these constraints making chargino responsible for large 
asymmetries. On the other hand, pure gluino exchanges can easily saturate both 
the constraints on B — > 4>K$ and B — > g'K$ CP asymmetries. Moreover, we also 
find that the simultaneous contributions from gluino and chargino exchanges could 
easily account for the present experimental results on the mentioned asymmetries. 
Remarkably, large experimentally allowed enhancements of B — > n'Ks branching 
ratio can easily be achieved by the contribution of two mass insertions in gluino ex- 
changes. Finally, we analyze the correlations between the CP asymmetries of these 
processes and the direct CP asymmetry in b — > S7 decay. When all experimental 
constraints are satisfied, supersymmetry favors large and positive values of b — ► 57 
asymmetry. 



1 Introduction 



The B-factories are producing interesting experimental results with continuously increas- 
ing integrated luminosities. Currently they offer one of the most promising routes to 
test the Kobayashi-Maskawa ansatz for the CP violation. The understanding of the CP 
violating mechanism is one of the major open problems in particle physics. There are 
reasons to believe that the Standard Model (SM) cannot provide a complete description 
for the CP violating phenomena in nature. For instance, it is established that the SM 
mechanism of CP violation cannot account for the observed size of the baryon asymmetry 
in the Universe, and additional sources of CP violation beyond the SM are needed [1]. 

Recently, BaBar and Belle collaborations [2] announced large deviations from the SM 
expectations in the CP asymmetry of B° — > 4>K S and branching ratio of B° — > rj'K . 
These discrepancies have been interpreted as possible consequences of new physics (NP) 
beyond the SM [3-12]. For B° and B° decays to a CP eigenstate fcp, the time dependent 
CP asymmetries are usually described by rates Q>f CP (t), 

a fcp (t) = r( fo (t) - fcp) - r(£? ° (t) - fcp) = C fcp cos AM B t + S fcp sin AM B t, (1) 

r(E°(f) - f CP ) + r(B°(t) -> f p) 

where Cf CP and Sf CP represent the parameters of direct and indirect CP violations re- 
spectively, and AMs d is the -B eigenstate mass difference. 

In the Standard Model, the angle (3 in the unitary triangle of Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [13], can be measured from B meson decays. The golden mode 
B° — > J/ipKs is dominated by tree contribution and Sj/^k s = sin 2/5 + 0(X 3 ), where A is 
the Cabibbo mixing angle (see e.g. [14]). These uncertainties are less than 1%. 

The dominant part of the decay amplitudes for B° — > <pK s , rj'Kg is assumed to come 
from the gluonic penguin [15, 16], but some contribution from the tree level b — > uus decay 
is expected. The \<f>) is almost pure \ss) and consequently this decay mode corresponds 
also accurately, up to terms of orders C(A 2 ), to sin 2(3 in the SM [17]. The b — > uus tree 
level contribution to B d — > rj'K was estimated in [16]. It was found that the tree level 
amplitude is less than 2% of the gluonic penguin amplitude. Thus also in this mode one 
measures the angle (3 with a good precision in the SM. Therefore, it is expected that NP 
contributions to the CP asymmetries in B° — > (f)Ks, rj' Kg decays are more significant 
than in B° — > J/ipKs and can compete with the SM one. 

New sample of data have been recently analyzed by BaBar and Belle collaboration 
and the results with higher statistics have been now announced in [18,19]. The new 
experimental value of the indirect CP asymmetry parameter for B° — > J/tpK s is given 
by [18, 19] 

Sj/^k s = 0.726 ± 0.037, (2) 

which does not differ much from the previous one [20,21], and agrees quite well with 
the SM prediction 0.715t[J;o45 [22]. However, results of Belle on the corresponding sin 2(3 
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extracted for B° — > (j)K s process has changed dramatically [18, 19] 

S 4)Ks = 0.50 ±0.25l&{2 (BaBar), 

= 0.06 ±0.33 ±0.09 (Belle), (3) 

where the first errors are statistical and the second systematic, showing now a better 
agreement than before [23,24]. However, as we can see from Eq.(3), the relative central 
values are still different. BaBar results [18] are more compatible with SM predictions, 
while Belle measurements [19] still show a deviation from the cc measurements of about 
2a. Moreover, the average S$k s — 0.34 ± 0.20 is quite different from the previous one 
[25], displaying now 1.7 'a deviation from Eq.(2). 

Furthermore, the most recent measured CP asymmetry in the B° — > rj'Ks decay is 
found by BaBar [18] and Belle [19] collaborations as 

S^Ks = 0.27 ±0.14 ±0.03 (BaBar) 

= 0.65 ±0.18 ±0.04 (Belle), (4) 

with an average S v >k s = 0.41 ± 0.11, which shows a 2.5cx discrepancy from Eq. (2). For 
the previous results see (BaBar) [26] and (Belle) [24]. 

It is interesting to note that the new results on s-penguin modes from both experiments 
differ from the value extracted from the cc mode (J/ip), BaBar by 2.7<r and Belle by 2.4<r 
[18,19]. At the same time the experiments agree with each other, and even the central 
values are quite close: 

0.42 ±0.10 BaBar, 0.43+°;}? Belle. 

On the other hand, the experimental measurements of the branching ratios of B° — > 
<pK° and B° -> rj'K at BaBar [27], Belle [28], and CLEO [29] lead to the following 
averaged results [25] : 

BR(B° -> <j)K ) = (8.3±i;g) x 10~ 6 , (5) 
BR(B° -> r)'K°) = (65.2±g;g) x lO" 6 . (6) 

From theoretical side, the SM predictions for BR(B — > <pK) 1 are in good agreement with 
Eq.(5), while showing a large discrepancy, being experimentally two to five times larger, 
for BR(B — > rj'K) in Eq.(6) [30]. This discrepancy is not new and it has created a 
growing interest in the subject. However, since it is observed only in B — > rj'K process, 
mechanisms based on the peculiar structure of meson, such as intrinsic charm [31] and 
gluonium [32] content, have been investigated to solve the puzzle. 

1 In order to simplify our notation, from now on everywhere, where the symbols B and K will appear, 
they will generically indicate neutral B° (B°) and K° (K°) mesons, respectively. 
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Supersymmetry (SUSY) is one of the most popular candidates for physics beyond the 
SM [33]. In SUSY models there are new sources of CP violation besides the CKM phase 
[34]. The soft SUSY breaking (SSB) terms contain several parameters that may be com- 
plex, as may also be the SUSY preserving \x parameter in the Higgs sector. Then, new- 
CP violating phases can naturally arise in the SSB sector of scalar- quarks (squarks) and 
-leptons (sleptons). These new phases have significant implications for the electric dipole 
moments (EDM) of the neutron, electron, and mercury atom [34, 35] and can be strongly 
constrained by the negative search of CP violation in EDM experiments. Therefore, it 
remains a challenge for SUSY to provide an explanation for the above mentioned dis- 
crepancies in I?-decays, while keeping the above EDMs within their experimental ranges. 
This can be possible, for instance, in SUSY models where the CP violating phases enter- 
ing in b — > s transitions are independent from the corresponding ones affecting EDMs. 
An interesting example is provided by SUSY models with flavor dependent CP violating 
phases [35-37]. 

In SUSY framework, the main effect on B — > <pK is usually assumed to come from 
the gluino loop contributions to s-penguin diagrams [6-9,38], if R parity is conserved. 2 
However, charginos could also be responsible for such discrepancy, as it has been discussed 
in [10, 11], although b — » 57 constraints strongly suppress their contribution [11]. Similarly 
to B — > (f>K, one would a priori expect large effects from SUSY corrections to B — > 
rj'K as well, which may contradict the experimental results reported in Eq.(4). Possible 
mechanisms to explain such behavior in the supersymmetric context have been proposed 
in Refs.[8,9,39]. 

In this paper we perform a complete analysis of SUSY contributions to the CP asym- 
metries and branching ratios of B — > <f)K and B — > rj'K processes. Previous analyses on 
the same issue have considered either gluino [6] or chargino exchanges [7, 10, 11]. We think 
that in the framework of general SUSY models a complete analysis involving both effects 
is in fact needed. Indeed, we find that a SUSY scenario in which both chargino and gluino 
give sizeable contributions to CP violating processes, is an interesting viable possibility. 
For instance, large effects of chargino contributions to CP asymmetries that one would 
expect to be excluded by b — > 57 constraints [11], could be achieved taking into account 
gluino exchanges. This is due to potentially destructive interferences between chargino 
and gluino amplitudes in b — > 57 decay that eventually relax b — > 57 bounds. 

In our analysis we consider both the effects of gluino and chargino exchanges by 
using the mass insertion method (MIA) [40]. As known, this method is a useful tool for 
analyzing SUSY contributions to flavor changing neutral current processes (FCNC) since 
it allows to parametrize, in a model independent way, the main sources of flavor violations 

2 Recently, SUSY effects to B — > <pK decays have been analyzed in the context of R-parity breaking 
models [5]. In this case, as well as in models with extra dimensions [4], effects are induced at tree-level. 
We will restrict our analysis to the R-parity conserving scenarios, where SUSY corrections to B — > <pK 
process always enter at 1-loop. 
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in general SUSY models. 

We take into account all the relevant operators involved in the effective Hamiltonian 
for AB = 1 transition, and provide analytical expression for the corresponding Wilson 
coefficients. We analyze the most interesting scenarios in which one or two mass insertions 
are dominant in both gluino and chargino sector. 

An important issue in these calculations is the method of evaluating the hadronic 
matrix elements for exclusive hadronic final states, which may play a crucial role in 
CP asymmetries. Many studies have been done with naive factorization approach (NF) 
[41, 42], for the computation of two-body nonleptonic B decays. Recently, a new approach 
has been developed, called QCD factorization (QCDF), [43-45], which offers the possi- 
bility to include non-factorizable contributions and to calculate the strong phases. The 
drawback in this approach is that it includes undetermined parameters ph,a and phases 
</>h,a, characterizing the infrared divergences. We critically consider both approaches and 
analyze theoretical uncertainties connected with SUSY predictions. We provide a com- 
parative study of SUSY contributions from chargino and gluino to B — > <pK and B — > rj'K 
processes in NF and QCDF approaches. We also analyze the branching ratios of these 
decays and investigate their correlations with CP asymmetries. Finally, we discuss the 
correlations between CP asymmetries of these processes and the direct CP asymmetry in 
b — > 57 decay [46]. 

This paper is organized as follows. In Section 2 we give the basic ingredients needed 
for calculations, with many important formulas provided in the Appendices. In Section 
3 we discuss the QCDF method in the evaluation of 5-meson decay amplitudes of B — > 
<pK and B — > rfK. Moreover, we compare SUSY predictions between NF and QCDF 
approaches. Section 4 is devoted to the CP asymmetries of these decay modes. Both 
effects from gluino and chargino exchanges are discussed. In section 5 the branching ratios 
and their correlations to asymmetries are considered. Section 6 contains the analysis of 
SUSY contributions to direct b — > S7 asymmetry, and correlations to the mentioned CP 
asymmetries are also considered. Finally, section 7 contains our main conclusions. 

2 SUSY contributions to the effective Hamiltonian 
of AB = 1 transitions 

We start our analysis by considering the supersymmetric effect in the non-leptonic AB = 1 
processes. Such an effect could be a probe for any testable SUSY implications in CP 
violating experiments. The most general effective Hamiltonian H^ B=1 for these processes 
can be expressed via the Operator Product Expansion (OPE) as [47] 
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+ {Qi^QiiQ^Ci} , (7) 

where X p = V p bV* s , with V p j, the unitary CKM matrix elements satisfying the unitarity 
triangle relation X t + A M + A c = 0, and Cj = Ci(fib) are the Wilson coefficients at low 
energy scale //& ~ 0{m b ). The basis = Qi(fib) is given by the relevant local operators 
renormalized at the same scale namely 

(Pabp)v~A {S/3P a )v~A 

{s a bp) V -A ^2(qpq a )v-A , 
q 

(s a bp)v-A ~y](Ql3Qa)v+A , 

q 

3 

(s a bp) V - A Y ^ e q(Wa)v+A , 
q A 

3 

= (s a bp)v-A Y ^ e q(.Wo)v~A , 
q 1 

^%V^(1+7 5 )GXV (8) 

Here a and (3 stand for color indices, and t£p the SU(3) C color matrices, = |i[7 M , 7"]. 
Moreover, e g are quark electric charges in unity of e, (??)y±A = 97/x(l ± 7s)?, an d Q runs 
over u, d, s, c, and b quark labels. In the SM only the first part of right hand side of 
Eq.(7) (inside first curly brackets) containing operators Qi will contribute, where Qi 2 
refer to the current-current operators, Q3-6 to the QCD penguin operators, and Q7-10 to 
the electroweak penguin operators, while Qj 7 and Q 8g are the magnetic and the chromo- 
magnetic dipole operators, respectively. In addition, operators Qi = Qi(fib) are obtained 
from Qi by the chirality exchange (q~iq2)v±A —> (qiQ2)v^a- Notice that in the SM the 
coefficients Ci identically vanish due to the V-A structure of charged weak currents, while 
in MSSM they can receive contributions from both chargino and gluino exchanges [48, 49]. 



QI = 


(Pb)v-A (sp) V -A , 


Q? = 


Qs = 


(sb) v - A J2(qq)v-A , 

q 


Qa = 


Q 5 = 


(sb) v ^ A J2(^)v+a , 
q 


Qe = 


Qi = 


3 

(sb)v- A Y ^e q (qq) v+A , 
q A 


Q 8 = 


Q9 = 


3 

{sb) V -A Y 7; e q(q<l)v-A , 
q 1 


Qio - 


Qii = 


■^m h sa" v {\ + 75) V, 


Q%g = 



Due to the asymptotic freedom of QCD, the calculation of hadronic weak decay ampli- 
tudes can be factorized by the product of long and short distance contributions. The first 
ones, that will be analyzed in the next section, are related to the evaluation of hadronic 
matrix elements of Qi and contain the main uncertainty of our predictions. On the other 
hand, the latter are contained in the Wilson coefficients Ci and they can be evaluated in 
perturbation theory with high precision [47]. For instance, all the relevant contributions 
of particle spectra above the W mass (m w ) scale, including SUSY particle exchanges, will 
enter in Ci(nw) at /iw ~ 0{m w ) scale. 

The low energy coefficients Ci(fib) can be extrapolated from the high energy ones 
Ci(nw) by solving the renormalization group equations for QCD and QED in the SM. 
The solution is generally expressed as follows [47] 

= Y Uijiv, Hw) Cj{nw), (9) 
3 
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where Uij(n, /^w) is the evolution matrix which takes into account the re- summation of the 
terms proportional to large logs (a s (/j, w ) log(/j, w / /j, b )) n (leading), (a 2 s (fi w ) \og(fi W / fj, b )) n 
(next-to- leading), etc., in QCD. 

In our analysis we include the next-to-leading order (NLO) corrections in QCD and 
QED for the Wilson coefficients C i= i-i as given in Ref. [47], while for C 77 (/i) and C$ g (fj,) 
we include only the leading order (LO) ones. 3 The expressions for the evolution matrix 
Uij(/i, at NLO in QCD and QED can be found in Ref. [47]. The reason for retaining 
only the LO accuracy in CV 7 (//) Cg g ([j,) is that the matrix elements of the dipole operators 
enter the decay amplitudes only at the NLO [43]. 

Next we discuss the SUSY contributions to the effective Hamiltonian in Eq.(7). The 
modifications caused by supersymmetry appear only in the boundary conditions of the 
Wilson coefficients at jiw scale and they can be computed through the appropriate match- 
ing with one-loop Feynman diagrams where Higgs, neutralino, gluino, and chargino are 
exchanged, see for instance Ref s. [48, 49, 51]. Only chargino and gluino contributions can 
provide a potential source of new CP violating phase in MSSM and could account for the 
observed large deviations in B — > <pKs asymmetry. In principle, the neutralino exchange 
diagrams involve the same mass insertions as the gluino ones, but they are strongly sup- 
pressed compared to the latter. For these reasons we neglect neutralino in our analysis. 
The charged Higgs contributions cannot generate any new source of CP violation in addi- 
tion to the SM one, or any sizeable effect to operators beyond the SM basis Qi. However, 
when charged Higgs contributions are taken into account together with chargino or gluino 
exchanges, their effect is relevant. In particular, as we will show in the next sections, due 
to destructive interferences with chargino and gluino amplitudes, b — > 57 constraints can 
be relaxed allowing sizeable contributions to the the CP asymmetries. 

The results for Wilson coefficients at /iw scale can be expressed as follows 

Ci(nw) = Cf + Cf + A7 1 {Cf + Cf} (10) 

where Cf , Cf 1 , Cf, and Cf correspond to the W, charged Higgs, chargino, and gluino 
exchanges respectively. In our analysis we will impose the boundary conditions for Cf , 
and Cf at the scale nw = m w, although they should apply to the energy scale at which 
SUSY particles are integrated out, namely M s . However, these threshold corrections, 
originating from the mismatch of energy scales, are numerically not significant since the 
running of a s from Ms to mw is not very steep. 

Finally, for NDR renormalization scheme, the electroweak contributions to the Wilson 
coefficients are given by [48, 51] 

w _ Ua s 

l07T 

3 However, in the evaluation of the BR of b — > decay, the complete NLO corrections in CV 7 (/Ltf,) have 
been taken into account [50] . 
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r (W,H, x ) a 1 ( B {W,H, X ) , 1 R (W,H, X ) , WW.ff.xA _ a « 

C (W,H, X ) = ^s_ E (W,H, x ) 
87T 

r (W,H, x ) a s p (W,H,x) 

5 ~ 24tt 

6 8tt 

C (W,H, X ) = ^_ ^ C (W,H, X ) + D (W,H, X )^ ; 

c 8 M = o! 

c (w,h, x ) = f AC {w,H, x ) + d (w,h, x ) + * f B Qv,H, x ) + B^ H ' x) - AC {w ' H ' x) S\ , 
67T V sin ' 

= 0, 

r r(w,/r,x) = m (w,h, x ) 

7 ' 

C (W,H, X ) = M (W,H, X ) (U) 

where a s and ct are evaluated at mw scale. The functions appearing above, include the 
contributions from photon-penguins (D), Z-penguins (C), gluon-penguins (E), boxes with 
external down quarks (Bd) and up-quarks (B u ), the magnetic- (M 7 ), and the chromo- 
magnetic-penguins (M g ). The corresponding SM results are [47] D w = D(x t ), C 1 ^ = 
C(x t ), -E 1 ^ = E(x t ), and B w = B(x t ) with = m^/rn^, and analogously for charged 
Higgs [51], £> H = Dh(x h ), C h = Ch(xh), £ h = E H (^), and 5" = B H (x H ) with 
Xij = mf/mjj, where the loop functions B, C, D, E and Bh, Ch, Dh, Eh and are provided 
in Appendix C. Regarding the SM and charged Higgs contributions to magnetic-penguins, 
we have [48] 

Mf = -x t (Fx(x t ) + j^te)) , Mf = -^F 1 (^) 
Mf = -^{^F 1 {x H )+F 2 {x H )^ C ot 2 p + ^F z {x H )+F A {x H ) 
Mf = -^-(F 1 (^)cot 2 /3 + F 3 (^' N 



where the functions F^x) are reported in Appendix C. Finally, the gluino and chargino 
exact contributions to the expressions appearing in Eq.(ll) can be found in Refs.[48, 51], 
while here we will provide only the corresponding results in the so called mass insertion 
approximation. 

Regarding the SUSY contributions to the opposite chirality operators Qt, we stress that 
while in the SM the Wilson coefficients C, identically vanish, in SUSY models chargino 
and gluino exchanges could sizeably affect these coefficients. However, in case of charginos, 
these effects are quite small being proportional to the Yukawa couplings of light quarks 
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[11] and so we will not include them in our analysis. Moreover, we have also neglected the 
small contributions to C* 2 coming from box diagrams, where both chargino and gluino 
are exchanged [11,51]. 

As mentioned in the introduction, in order to perform a model independent analysis of 
FCNC processes in general SUSY models, it is very convenient to adopt the mass insertion 
approximation (MIA) method [40]. This method has been applied in many analyses of 
FCNC processes in K, D, and B meson sector and leptonic sector, mediated by gluino 
and neutralino exchanges, respectively [40]. More recently, it has been extended to the 
sector of FCNC processes mediated by chargino exchanges [11,52,53]. In particular, in 
[11], the chargino functions appearing in Eq.(ll) have been calculated at the first order 
in MIA. 

Let us briefly recall the main steps of this approximation. In MIA framework, one 
chooses a basis (called super-CKM basis) where the couplings of fermions and sfermions 
to neutral gaugino fields are flavor diagonal. In this basis, the interacting Lagrangian 
involving charginos is given by 

W = -9 EE {Vki KL dl (xtT u b L - U* k2 {Yt g -K + ) ab d» (xtT ul 

k a,b 

-V k2 (K.Y u d ^) ab d a L (xtTu b R ), (12) 

where q^i — |(1 ± 75)?, and contraction of color and Dirac indices is understood. Here 
Y tt d ^ ag are the diagonal Yukawa matrices, and K stands for the CKM matrix. The indices 
a, b and k label flavor and chargino mass eigenstates, respectively, and V, U are the 
chargino mixing matrices defined by 

w - •»*>. - ( ^T n ' ) • (13) 

where M 2 is the weak gaugino mass, ji is the supersymmetric Higgs mixing term, and tan (3 
is the ratio of the vacuum expectation value (VEV) of the up-type Higgs to the VEV of 
the down-type Higgs 4 . As one can see from Eq.(12), the higgsino couplings are suppressed 
by Yukawas of the light quarks, and therefore they are negligible, except for the stop- 
bottom interaction which is directly enhanced by the top Yukawa (Y t ). In our analysis we 
neglect the higgsino contributions proportional to the Yukawa couplings of light quarks 
with the exception of the bottom Yukawa Y b , since its effect could be enhanced by large 
tan (3. However, it is easy to show that this vertex cannot affect dimension six operators of 
the effective Hamiltonian for AB = 1 transitions (operators <5i=i-io in Eq.(7)) and only 
interactions involving left down quarks will contribute. On the contrary, contributions 
proportional to bottom Yukawa Y b enter in the Wilson coefficients of dipole operators 
(CV 7 , Cs g ) due to the chirality flip of b — > 57 and b — > sg transitions. 

4 This tan/3 should not be confused with the angle (3 of the unitarity triangle. 
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Figure 1: Penguin diagrams for AS = 1 transitions with chargino (xj) exchanges at the first 
order in mass insertion. Here U, U' = {u, c, t}, with indices h, k,m,n = {L, R} and I, J = {1, 2}. 
The cross symbol in the squark propagator indicates the mass insertion. The corresponding 
diagrams at zero order in mass insertion are simply obtained by removing the mass insertion in 
the propagators of up- type squarks (U). 
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Figure 2: Box diagrams for AB = 1 transitions with chargino exchanges at the first order in 
mass insertion, where U,U' = {t,c,u}, U = {c, u}, and D = {b, s,d}, where h,k,l,n,r,m,q = 
{L,R}. 
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Figure 4: Box diagrams for AB = 1 transitions with gluino exchanges at the first order in mass 
insertion, where q = {b, c, s, d, u} and h,k,m = {L, R}. 
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As mentioned in the case of MIA, the flavor mixing is displayed by the non-diagonal 
entries of the sfermion mass matrices. Denoting by A 9 AB the off-diagonal terms in the 
sfermion (q = u, d) mass matrices for the up and down, respectively, where A, B indicate 
chirality couplings to fermions A,B — (L, R), the A-B squark propagator can be expanded 

as 

< a «S> = i (fl - *»1 - AS,)"' = j*^ + J^U + 0(A2) , (14) 

where q = u,d selects up or down sector, respectively, a,b = (1,2,3) are flavor indices, 
1 is the unit matrix, and m is the average squark mass. As we will see in the following, 
it is convenient to parametrize this expansion in terms of the dimensionless quantity 
(S q AB ) a b = (A q AB ) ab /m 2 . At the first order in MIA, the penguin and box diagrams which 
contribute to the AB = 1 effective Hamiltonian are given in Figs. 1 and 2, respectively. 
Evaluating the diagrams in Figs. 1 and 2 by retaining only terms proportional to bottom- 
and top-quark Yukawa couplings and performing the matching, the chargino contributions 
to the Wilson coefficients in Eqs.(ll) can be determined from the following relations [11] 



a,b a 



J RL)3a 



RL 



Y t R F 



+ 



[ £ Kl 2 K aZ {5l R ) a ,] Y t R L F R + [K* 2 K 33 ((8 u rr ) 33 R rr + R° F ) ] Y t 2 , (15) 



where the symbol F = {D, C, E, B d: B U) M 7 , M g }, and the detailed expressions for R F L , 
R F R , R F L , R F R , and R F can be found in Appendix A 5 . Notice that the last term in 
Eq.(15), proportional to R F , is independent of mass insertions. This is due to the fact 
that for chargino exchanges the super-GIM mechanism is only partially effective, when 
only squarks but not quarks are taken to be degenerate. 

Here we will just concentrate on the dominant contributions which turn out to be 
due to the chromo-magnetic (M 9 ) penguin and Z-penguin (C) diagrams [11]. From the 
above expressions, it is clear that LR and RR contributions are suppressed by order A 2 
or A 3 , where A = sin# c ~ 0.22, with 9 C the Cabibbo angle. In our analysis we adopt 
the approximation of retaining only terms proportional to order A = sin 8 C . In this case, 
Eq.(15) simplifies as follows [11] 

F x = ^ LL R L F L + Y t ^ RL Rf', (16) 

where £ LL = (5't L ) 32 + A {8 u LL ) 3l and £ RL = {5 u RL ) m + A (5% L ) 3V 

The functions R F L and R F L depend on the SUSY parameters through the chargino 
masses (m Xi ), squark masses (m) and the entries of the chargino mass matrix. For instance 
for Z and magnetic (chromo-magnetic) dipole penguins Rq L ' rl and Rm'^, respectively, 

5 The expression R F was missing in [11]. The contribution is 0{\ 2 ) and does not change the numerical 
results. 
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we have 



R L c L = £l^| 2 4 0) (^)+ E [UaVnU^pW^Xj 

i=l,2 i,j"=l,2 



+ |Vy 2 |^| 2 Q-P^,^ 



^ = 4E K 2 ^i4 0) (^)- E ^(^P^,^) 

Z i=l,2 i,j=l,2 
+ V^ 1 P£ ) (^)), 



= T,\^i\ 2 ^P^ g (x t )-Y b Y / V ll U l2 x w ^P^ g (x l ), 

i i m b 

R mL = -EVi^xwiPEttix,), (17) 

i 

where x W i = rn^/m^., Xi = rri^./m 2 , Xi = m 2 /m 2 ., and Xij = rri^./m 2 ,.. The loop 

functions Pc' 2 \x,y), P^^ R \x) are provided in Appendix C. Finally, U and V are the 
matrices that diagonalize the chargino mass matrix, defined in Eq. (13). 

It is worth mentioning that the large effects of chargino contributions to Cy 7 and C 8g 
come from the terms in P^ and R^f , respectively, which are enhanced by m Xi /mb in 
Eq.(17). However, these terms are also multiplied by the Yukawa bottom Yb, which leads 
to enhancing the coefficients of the LL mass insertion in CV 7 and C§ 9 at large tan/3. 
As we will see later on, this effect will play a crucial role in chargino contributions to 
B — > (p(r]')K decays at large tan/5. 

We will also consider the case in which the mass of stop-right (irii R ) is lighter than 
other squarks. In this case the functional form of the Rp remains unchanged, while only 
the expressions of Rp L should be modified in the way described in Appendix B. 

Now let us turn to the gluino contributions in the b — > s transition. In the super-CKM 
basis, the quark-squark-gluino interaction is given by: 

C dtg = V2g s T A {(d a P L ~g A ) d R - (rf» P R g A ) d{ + h.c] , (18) 

where g A are the gluino Majorana fields, dr RL the squark fields, T A are the SU (3) c gen- 
erators, and a, (3 are color indices. The dominant gluino contributions are due to the 
QCD penguin diagrams, and the magnetic and chromo-magnetic dipole operators. At the 
first order in MIA, the penguin and box diagrams are shown in Figs. 3 and 4, respec- 
tively. Performing the matching, the gluino contributions to the corresponding Wilson 
coefficients at SUSY scale are given by [49] 6 



'2V2G F mj} ULLj ™ 



gBl (x) - -B 2 (x) - -P^x) - -P 2 (x) 



6 Note that the Wilson coefficients of Eq.(19) are different from those reported in Ref. [49] by a minus 
sign and a rescaling factor. This is due to the different convention for the Wilson coefficients in the 
effective Hamiltonian of Eq.(7). 
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5 — 



Cf 



6 — 



o ?7 



r<9 



9V2G F mj 



2V2G F ml K LL) ™ [ 3 v ; 3 v ; 6 v ; 2 A 
— ^ f /) f — Si fx) + — B 2 (x) - — Pi (x) - -P 2 

M 3 (i) + 3M 4 ( I )) + (4 



V2G F m? 



d 

LL 



1 

23 V 3' 



,(19) 



where Ci,8g are obtained from C^gg by exchanging L <-> P in (5^)23. The functions 
appearing in these expressions can be found in appendix C, with a; = m~/ra~. 

Now we would like to comment about the chiral enhancement in Cj Jj8g due to the 
(J>lr) 23 mass insertion. As for chargino in Eq.(17), the term proportional to (J>lr) 



23 



in Cj jSg in Eq.(19), has also the large enhancement factor rrig/mb in front. Moreover, 
contrary to the chargino case, this term is not suppressed by the bottom Yukawa coupling. 
As we will see later on, this enhancement factor will be responsible for the dominant gluino 
effects in B — > 4>(rj')K decays. 

In concluding this section, we emphasize that in the SM the b — > s transition process 
is dominated by the top quark mediated penguin diagram, which does not include any 
CP violating phase. Therefore, this process, as observed at the B-factories, opens up the 
possibility to probe virtual effects from new sources of flavor structure and CP violation. 
The SUSY contributions through gluino and chargino exchanges are independent. The 
ones from gluino depend on the flavor structure of the down squark sector, namely (6^)23, 
while the other ones from chargino depend on the up squark sector, particularly (6^3)32 
and (5^5)31. So, depending on the constraints imposed on the flavor structure of the 
down or up sector (for instance from b — > 57 decay, K — K and B — B mixing), gluino 
or chargino exchanges could give sizeable effects. As known, in many SUSY scenarios the 
lighter chargino is expected to be one of the lightest supersymmetric particles. Thus, it 
could contribute significantly in the one-loop processes. However, even though gluino in 
most models is expected to be heavier than chargino, it is a strongly interacting particle, 
and may give the dominant effect as well. 



3 B —> (f>(r} f )K in QCD factorization approach 

The calculation of B — > 4>(r]')K decays involves the evaluation of the hadronic matrix 
elements of related operators in the effective Hamiltonian, which is the most uncertain 
part of this calculation. In the limit in which m b ^> Aq C d and neglecting QCD corrections 
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in a s , the hadronic matrix elements of B meson decays in two mesons can be factorized, 
for example for B — > MiM 2 , in the form 

(M 1 M 2 \Q t \B°) = (M^B ) x (M 2 \j 2 \0) (20) 

where M 12 indicates two generic mesons, Qi are local four fermion operators of the ef- 
fective Hamiltonian in Eq.(7), and ji i2 represent bilinear quark currents. Then, the final 
results can be usually parametrized by the product of the decay constants and the tran- 
sition form factors. This approach is known as naive factorization (NF) [41,42]. Then, 
the hadronic matrix elements for B — > <pK are given by [42], 

(^°|Qi|S°) = 0, (<PK°\Q 2 \B°) = 0, (<PK°\Q 3 \B°) = |x, 

(0^°|g 4 | J B ) = ((f>K \Q 5 \B ) = X, (<PK°\Q 6 \B°) = ^X, 

{4>K Q \Q 7 \B G ) = \e a X, (<f ) K°\Qs\B°) = l-e s X, (<pK°\Q 9 \B ) = 3e s X, 



2 ^ ^ 2 

(0A-°|gio|5°) = e s X, (0#°|Q 77 |fl°) = 0, 



(21) 



where for color number N c = 3 one gets 

X = 2F^ K (ml)f <p m^p K • e ). (22) 

In Eq.(22), is the meson mass, F^~" K (m^) is the transition form factor evaluated 
at transfered momentum of the order of scale, Pk stands for K momentum, and 
is the polarization vector. In our analysis, for the parameters above, we will use the 
central values = 1.02 GeV, / = 0.233 GeV, and F^ K = 0.35 ± 0.05, and for the 



scalar product (p K ■ = ^ [_JL_ ( m | - + mj)] - ~ 13 GeV. 

Since the evaluation of the matrix element of Q$ g goes beyond the standard NF ap- 
proach, we report here its result as given in [42] 



(<f>K \Q8 9 \B°) - " ""' 



(Qa) + (Q 6 ) - I «Q 3 > + W) 



(23) 



where (Qi) = (0X°|<5i|-B°). In the derivation of Eq.(23), the following assumption was 
made [42] 

q* 1 = v^— • (24) 

The momentum g M appearing in Eq.(24), is connected to the virtual gluon in the Qs g 
operator, and it is given by g M = p^—p^-, with the momentum of b(s) quark, and {q 2 ) 
is an averaged value of q 2 . It has been shown that the physical range of q 2 in B — > 
is m 2 /4 < (g 2 ) < m 2 /2 [54]. Notice that the term 1/ J (q 2 ) in Eq.(23) has origin from the 



propagator of the virtual gluon exchange between Q 8g and external sources. 
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Here we stress that in the SM, the large uncertainty on (q 2 ) does not necessarily con- 
vert in a large uncertainty on the B — > <pK amplitude, since the contribution of Q$ g is 
not the dominant source in the SM. On the contrary, (q 2 ) could cause a large uncertainty 
in the numerical analysis of SUSY models. Indeed, as we will show in the next section, 
in most relevant SUSY scenarios, C 8g provides the dominant source to B — > <pK ampli- 
tude. Moreover, the NF approach suffices a serious problem, namely, the decay amplitude 
in this approximation is not scale independent. The hadronic matrix elements cannot 
compensate for the scale dependence of the Wilson coefficients. This might be a hint for 
the necessity of including higher order QCD corrections to the hadronic matrix elements. 
Furthermore, due to the approximations used in NF, one cannot predict the direct CP 
asymmetries due to the assumption of no strong re-scattering in the final state, thus leav- 
ing undetermined the predictions of strong phases. 

Recently, in the framework of QCD and heavy quark effective theory, a more consis- 
tent method for the determination of nonleptonic B meson decays has been developed 
[43]. In this approach, called QCD factorization (QCDF), the hadronic matrix elements 
can be computed from the first principles by means of perturbative QCD and hQCD/ m b 
expansions. Final results can be simply expressed in terms of form factors and meson 
light-cone distribution amplitudes. Then, the usual NF is recovered only in the limit in 
which AQ C D/ m b and a s corrections are set to zero. A nice feature of QCDF is that the 
strong phases of non-leptonic two body decays can be predicted. 

In QCDF the hadronic matrix element for B — > M K with M = 0, rj in the heavy 
quark limit m;, ^> Aqcd can be written as [43] 

(MK\Q l \B) QCDF =(MK\Q l \B) 

NF- 

where {MK\Qi\B) NF denotes the NF results. The second and third term in the bracket 
represent the radiative corrections in as and Agc^/m^. Notice that, even though at 
higher order in a s the simple factorization is broken, these corrections can be calculated 
systematically in terms of short-distance coefficients and meson light-cone distribution 
functions. 

Now we briefly recall the main results of this method [43,44]. In QCDF the decay 
amplitudes of B — > (p(rj')K can be expressed as 

A (B -> <j>{rf)K) =A f (B^ 4>{r]')K) + A a (B -> <j>{rf)K) , (26) 

where 

ft 10 

Af (B -> </>(rf)K) = ^EE V pb V p * s at {T,,) (<f>(ri')K\Qi\B)NF, (27) 

V ^ p=u,c 1=1 

and 

A a (b - Mm = %fBf K n E E y Pb y; s t v,) - m 

p=u,c i=l 



n. \ 



A 



QCD 



(25) 
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The first term *4/ (B — > (f>{rf)K) includes vertex corrections, penguin corrections and hard 
spectator scattering contributions which are involved in the parameters af^ v \ The other 
term A a (B — > (j>(rf')K) includes the weak annihilation contributions which are absorbed 
in the parameters bf^\ However, these contributions contain infrared divergences, and 
the subtractions of these divergences are usually parametrized as [44] 

f ld -E^X HiA = (l +PHiA e^yn(^-), (29) 

Jo X V 1 \AqcdJ 

where Ph,a are free parameters expected to be of order of Ph,a — and </>h,a £ [0, 27r]. 

As already discussed in Ref.[44], if one does not require fine tuning of the annihilation 
phase (f>A, the pa parameter gets an upper bound from measurements on branching ratios, 
which is of order of pa < 2. Clearly, large values of pa are still possible, but in this 
case strong fine tuning in the phase (pA is required. However, assumptions of very large 
values of Ph,a, which implicitly means large contributions from hard scattering and weak 
annihilation diagrams, seem to be quite unrealistic. 

Notice that the annihilation topology contribution due to the effective operator Q 8g 
has not yet been calculated, and it is expected to have as well a logarithmic divergence. 
This effect should increase the theoretical uncertainty substantially specially in models 
like supersymmetric ones, where the Qs g plays a crucial rule in the B — > (f)K decay. 

Following the scheme and the notation of Ref. [44] we write the decay amplitude of 
B — > <pK as: 

A(B -> (j>K) = -iA K;4 , X p [ a 3 + "4 - \°%,ew ~ \ a l,EW + 03 

p=u,c 

~ 2^1,ew + Ps3 ~ 2^S3,EW 5 (30) 

where \ p = V p bV* s and A K ^ = ^m 2 B f<f,F+~* K . The quantities af = af (K, 0) and 
f3f = (3f(K, 0) depend 7 , in addition to the leading contribution of the b — > s transition, 
on the one loop vertex corrections, hard spectator interactions and penguin corrections 
[44]. 

As mentioned in section 2, NP effects are parametrized in the Wilson coefficients, 
while all the other functions involved in the definition of af and (3f depend on some 
theoretical input parameters like the QCD scale Aqcd, the value of the running masses, 
and parameters of vector-meson distribution amplitudes. Therefore it would be very useful 
for future analyses involving any NP scenarios, to provide a numerical parametrization 
of Eq.(30) in terms of the Wilson coefficients at low energy. In this respect it is very 
convenient to define new Wilson coefficients Cj and C, according to the parametrization 

7 Notice that in the notation of Ref.[44], the order of the arguments in the functions af (Mi,M<2>) and 
/3f (Mi, M2) is fixed and it is determined by the order of the arguments (Mi, M2) in the pre-factor Am l ,m 2 , 
where Mi >2 labels the final states. In Eq.(30) this order corresponds to Mi = K, M 2 = 4>. 
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of the effective Hamiltonian in Eq.(7) 



G 



(31) 



where the operators basis Q { and Qi are the same ones of Eq.(8). 8 Fixing the experimental 
and SM parameters to their center values as given in Table 1 of Ref.[44], we can present 
the explicit dependence of the decay amplitude on the Wilson coefficients relevant for this 
process. In particular, for 

A(B -> 4>K) = -i^m 2 B F^ K U ]T #(0)(Q + Q), (32) 



we obtain 



i=l. .10,77,! 



# 4 (0 
#5(0 

# 6 (0 

H 7 ((f> 

# 8 (0 
#9(0 
#10(0 



H- 



7j{ 



-0.0002 - 0.0002i 
0.011 + 0.009i, 

-1.23 + 0.089i - 0.005X A - O.OOOGX^ _ qqi 3Xh , 

-1.17 + 0.13i-0.014X H , 

-1.03 + 0.053i + 0.086X A - 0.008X1 

-0.29 - 0.022i + 0.028X A - 0.024X A + 0.014X^, 

0.52 - 0.026i - 0.006X A + 0.004X A , 

0.18 + 0.037i - 0.019X A + 0.012X A - 0.007XH, 

0.62 - 0.037i + 0.003X A + 0.0003X A + 0.007X H , 

0.62 - 0.037i + 0.007X H , 

-0.0004, 

0.047. 



(33) 



Notice that here both H 7j and H 8g do not have any dependence on X# A , since, as we 
mentioned, the hard scattering and weak annihilation contributions to Q 7l and Q 8g have 
been ignored. 

As can be seen from expressions in Eq.(33), terms proportional to Xa,h represent 
always small corrections for typical values of X# A ~ 0(1). Moreover, if we set X A # to 
zero, the contribution to the amplitude is (incidentally) quite close to the corresponding 
one in NF for (q 2 ) = m\j '4. Remarkably, the strong phases appearing in the terms 
independent on Xh,a in -£/j(0) are also negligible. Thus, in the limit Xh,a —> 0, the NF 
result (where strong phases are assumed to vanish) seems to be recovered for the choice 
(q 2 ) = ml/ 'A. 

8 Notice that in the notation of [44], the operators Qi and Q 2 correspond to our operators Q2 and Qi 
respectively. 
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Note that largest corrections in Xa are contained in the H 6 (<j)) term. In particular, for 
values of pa = l and 4>a — 0, we have Xa — 4 and so the term proportional to X\ in H e (cj)) 
becomes of the same order as the other term, which is independent on Xa- Therefore, 
it is expected that the effect of the weak annihilation parameter would be important 
when contributions to C e becomes large. Clearly, this last consideration applies in a NP 
context, where the new contribution to Wilson coefficients Cj could sizeably differ from 
the corresponding SM ones. 

Finally we would like to comment about the fact that contributions from Cj and Cj 
to the decay amplitude A(B — > <pK) are identically the same (with the same sign). This 
can be simply understood by noticing that 

(4>fP\Qi\B°) = (^Q^B ) . (34) 

which is due to the invariance of strong interactions under parity transformations, and to 
the fact that initial and final states have same parity. Indeed only terms proportional to 
the chiral structure V x V or A x A in the four fermion operators of Qi basis contribute 
to the matrix elements, and consequently property in Eq.(34) easily follows. Analogous 
considerations apply to Q 7l and Q 8g as well. 



Now we turn to the B — > rj'K process. As known, the physical state of rj and 7/ are 
the mixture of SU(3) singlet r] s = \ss) and octet r] q = (\uu) + \dd)^j components: 



\r]) \ _ f cos -sin0 \ / \r) q ) 
|V) / I sin0 cos0 J \ 1 77 s ) 



(35) 



Therefore, we have the following decay constants: f^, = f q sin <fi and ft = f s sin <fi and 
from a fit to experimental data one finds [44]: /, = (1.07 ± 0.02)/ w , f s = (1.34 ± 0.06)/ w 
and = 39.3 ±1.0. 

In the NF approach, the hadronic matrix elements of the B — > rj'K process are given 
by [9,42] 



(V'K°\Q 3 
WK°\Q 4 
(V'K°\Q 5 



(r/K°\Qt 



WK°\Q 7 



-X 2 (r ] 'K°\Q 2 \B°}=X 2 

l -X, + 2X 2 + ^X 3 

2 4 
Xi -\- — X 2 -\- —X 3 

3 3 

?±X X - 2X 2 - (l - 
B°) = - jjx 2 - 



B°) 
B°) 
B°) 
B°) 



B°) 



3 



X, 



- - R2 )X: 



R1X1 



-x 2 +ii--z\x. 
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{rfK°\Q s \S°) 
{r, 1 K a \Q 9 \B°) 
(r)'K°\Q 10 \B ) 



with 



X 1 

x 2 
x 3 
Ri 



-^ + x 2 - A -x 3 
-x 1 + ^- A -x 3 

3 3 



Xw 



-Ri\X: 



(36) 



<rn%-m*)F^(m 2 K )-jLf K , 



\m, 



m 



2m 2 K 



«>Fr K H,)pp K -p n Y v ,, 

{2m\ - ml) 



(m b - m s )(m s + m d ) ' 



Ri = 



(m b - m s )m s 



where F x (g ) = 0.35 is the B — n transition form factor evaluated at q scale, and 
Ik(-k) = 0.16(0.13) GeV is the decay constant of K(ir) meson. X v ' = 0.57 and = 0.82, 
which correspond to 9, p = —20°, represent the rate of the uu + dd and ss component in the 
rf . These matrix elements show that the B — > rj'K process receives a small contribution 
from color suppressed b — > uus tree diagram, in addition to the b — > sqq (q = u, d, s) 
penguin diagrams. As in case of B — > (pK decay, for the matrix element of the chromo- 
magnetic operator we have 



(rj'K^Q^B ) = 



a s m b 
47T ^/q 1 



(Qa) + (g 6 > - 3 «Q 3 ) + (<55» 



(37) 



In the QCDF approach, the decay amplitude of B — > rj'K is given by [44, 45] 



A(B^r)'K) = -^=A K , v , q {X P [S pu a 2 + 2a p 3 + -(^,ew + 

p=u,c " 



2Ps3 ~ Ps3,EW 



+ iA K , v > s [«3 + «4 - y a 3,EW ~ 7j a l,EW + 03 ~ ^3,EW + @S3 ~ yfls3,EW 



+ iA K ^ c 5 pc a 2 + «3 + -^ A v' q ,K 



>,p _ l c ,p 
2 



a 4 o a 4,EW 



+ 



P P 3-l(3 P 3,Ew}}, (38) 



where 9 



A 



K „, = ^£mlF B - K r v , ~ 6.8 xKT 6 GeV, 



9 Same notation as in A{B — > has been adopted here for and /3f, where = cuf (Mi, M2), 
/3f = <(Mi,M 2 ). The expressions for of (A", 77'), «?(»?', -^), P P i{K,i), and (3f(ri',K) can be found in 
Ref.[44,45]. 
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A 



K,rf c 



71 B J "' 

GF -mlF B ^ 



A n , 



V2 

G p 



1.05 x 10~ 5 GeV, 
-2.3 x 10~ 7 GeV, 



^mlF B ^'f K ~ 7.0xl0~ 6 GeV, 



(39) 



where the third contribution is negative due to the definition of negative decay constant 
f£i in [45] . As in the case of B — > <f)K, we will provide below the numerical parametrization 
in terms of the Wilson coefficients defined in Eq.(31). By fixing the hadronic parameters 
with their center values as in Table 1 of Ref.[44], we obtain 

.G F 
■ — ='* 
y/2 



A(B 



rfK) = -i^=m 2 B F^ K 



frf 



E 

=1.. 10,77,,' 



Hi(v')(Ci - Ci), 



(40) 



where 







0.44 + 0.0005i, 


H2W) 


r^j 


0.076 - 0.064i + 0.006X H , 




r^j 


2.23 - 0.15i + 0.009X A + 0.0008X1 + 0.014X H , 




r^j 


1.76- 0.29^ + 0.026^, 


H B (rf) 


r^j 


-1.52 + 0.004X A + 0.008X1 


HeW) 


r^> 


0.54 - 0.29i + 0.006X A + 0.027X1 + 0.026X H , 


HrW) 


r^j 


0.078 + 0.001X A - 0.004X1, 


HsW) 


r^j 


-0.58 + 0.02i + 0.004X A - 0.014X1 - 0.004X^, 


HM) 


r^j 


-0.44 + 0.054i + 0.005X A - 0.0004X1 - 0.007X 


Hio(v') 


r^j 


-0.80 + 0.02i - 0.004X H , 


H 7y (rf) 


r^j 


0.0007, 


H 8gW) 




-0.089. 



(41) 

The sign difference between Q and Cj appearing in Eq.(40) is due to the fact that, 
contrary to the B — > <pK transition, initial and final states have here opposite parity. 
Thus, due to the invariance of strong interactions under parity transformations, only 
V x A or A x V structures of four-fermion operators will contribute to the hadronic 
matrix elements, and so 

(r t 'K\Q i \B) = -(r t 'K\Q i \B). (42) 

The comparison between the coefficients Hi((f>) and H^r]') in Eqs.(41) tells us that, 
apart from the sign difference of the Wilson coefficients Q in the amplitudes of B — > rj'K 
and B — > <j)K, in QCDF these amplitudes are expected to be different, unlike the case of 
NF as discussed in Ref.[9]. However, notice that the values of Hi((f>) and Hi{rj) are quite 
sensitive to the scale \i where they have been evaluated. Our results in Eqs.(33),(41) 
correspond to the choice \x = m£ olc = 4.5 GeV. 
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4 CP asymmetry of B — > (f>(r] f )Ks: gluino / chargino 

Here we analyze the supersymmetric contributions to the time dependent CP asymmetries 
in B — > 07^5 and £> — > 77' Kg decays in the framework of mass insertion approximation, 
in gluino and chargino dominated scenarios. 

New physics could in principle affect the B meson decay by means of a new source of CP 
violating phase in the corresponding amplitude. In general this phase is different from the 
corresponding SM one. If so, then deviations on CP asymmetries from SM expectations 
can be sizeable, depending on the relative magnitude of SM and NP amplitudes. For 
instance, in the SM the B — > <pK s decay amplitude is generated at one loop and therefore 
it is very sensitive to NP contributions. In this respect, SUSY models with non minimal 
flavor structure and new CP violating phases in the squark mass matrices, can easily 
generate large deviations in the B — > <pKs asymmetry. 

As mentioned in the introduction, the time dependent CP asymmetry for B — > <pKs 
can be described by 

T(B (t) -> <pK s ) + T(B(t) -> (pK s ) 

(43) 

where C^k s and S^Ks represent the direct and the mixing CP asymmetry, respectively 
and they are given by 



_ \p(<f>K s )\ 2 -l Q 2Im[lp{4>K s ) 



\p{<t>K s )\* + r \p{<t>K s )\* + i 

The parameter p((pK s ) is defined by 



(44) 



where A((j)Ks) and A{(j>Ks) are the decay amplitudes of B° and B° mesons, respectively. 



Here, the mixing parameters p and q are defined by \Bi) = p\B) + q\B ), \B 2 ) = 
p\B) — q\lf) where |-Bi( 2 )) are mass eigenstates of B meson. The ratio of the mixing 



1 .-xe.Vjyu 



parameters is given by 

1 = -e-""^, (46) 

V v tb v t d 

where 9d represent any SUSY contribution to the B — B° mixing angle. Finally, the above 
amplitudes can be written in terms of the matrix element of the AB = 1 transition as 

A^Ks) = (<PK s \H* B f = ^E ), A^Ks) = (<t>K s \ (H$=*f \B°). (47) 
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In order to simplify our analysis, it is useful to parametrize the SUSY effects by 
introducing the ratio of SM and SUSY amplitudes as follows 

/ 4 SUSY \ 

— = R+e»*e«* (48) 

\ A / 4>K S 

(49) 

and analogously for the rj'Ks decay mode 

/ A SUSY \ 

KiH - Rte^e"* (50) 

where Ri stands for the corresponding absolute values of | a Jsm I ; the angles 9^ ^ are 
the corresponding SUSY CP violating phase, and 5^ ^ = S^^L — ^ U ^, Y parametrize the 
strong (CP conserving) phases. In this case, the mixing CP asymmetry S^k s m Eq.(43) 
takes the following form 

sin 2(3 + 2R 4> cos ^ sin(^ + 2(3) + R\ sin(2^ + 2(3) 
^ Ks l + 2R 4> cosS 4> cos6^ + Rl ' { } 

and analogously for B — > i]'K s 



sin 2(3 + 2R V , cos 5 V > sin(^ + 2(3) + R\, sin(2^ + 2(3) 
v ' Ks 1 + 2R V , cos 5^/ cos 9^ + ' ^ ^ 

Assuming that the SUSY contribution to the amplitude is smaller than the SM one 
i.e. R(j>(rf) ^ 1) one can simplify the above expressions as: 

S<t>(n')K s = sin 2/9 + 2 cos 2(3 sin 6^/) cos S^R^) + 0(-R^ (7/) ) . (53) 

It is now clear that in order to reduce S^Ks smaller than sin 2(3, the relative sign of sin.9^ 
and cos<50 has to be negative. If one assumes that sin 9^ cos 5^ ~ —1, then R^ > 0.1 is 
required in order to get S<f,K s within la of the experimental range. 

We will work in a minimal scheme of CP and flavor violation. This means that, in 
the framework of MIA, we will assume a dominant effect due to one single mass insertion. 
In this case the CP violating SUSY phase will coincide with the argument (module ir) of 
corresponding mass insertion. Moreover, in the following we will generalize this scheme 
by including two (dominant) mass insertions simultaneously, but assuming that their CP 
violating phases are the same. We will perform this analysis in both gluino and chargino 
scenarios. 10 

In the following subsections we present and discuss our results separately for B — > 4>Ks 
and B — > r]'K s processes. 

10 There are corrections to the effective Hamiltonian AB = 1, mediated by both chargino and gluino 
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4.1 CP asymmetry in B — > (j)Ks 

As shown in Eq.(3), the most recent world average measurement of the CP asymmetry 
S^k S i indicates a 1.7 a deviation from the cc measurement. In particular, 

S 4>Ks = 0.34 ± 0.20 (54) 

allows also for negative values of CP asymmetry at 2a level. In the light of these results, 
we will analyze here, in a model independent way, the SUSY scenarios, which are more 
favored to explain such deviations. Let us start first with gluino contributions to the CP 
asymmetry in B — > <pKs- 

These effects have been analyzed in Refs. [6, 8, 9] in the framework of NF or QCDF and 
by adopting the MIA method. However, from these works a direct comparison between 
the results in NF and QCDF predictions is not always possible, mainly because of different 
approaches used in analyzing the SUSY contributions. For this reason we reanalyze here 
the gluino contributions in both NF and QCDF approaches in a unique SUSY framework, 
by comparing the predictions for the CP asymmetries versus the relevant SUSY CP vi- 
olating phases. In QCDF we have included the last updated results for hard spectator 
scattering and weak annihilation diagrams [44], as explained in section 3. 

We present our numerical results for the gluino contributions to CP asymmetry S$k s 
in Figs. 5. The left and right plots correspond to the evaluation of amplitudes by means 
of NF and QCDF methods, respectively. In all the plots, regions inside the horizontal 
lines indicate the allowed 2a experimental range. In the top and bottom plots only one 
mass insertion per time is taken active, in particular this means that we scanned over 
I < 1 an d \ ($lr) I < 1- Then, SW S is plotted versus 0$, which in the case of 

V / 23 » * 23 

one dominant mass insertion should be identified here as 6^ = arg[(<5^ B )y]. 

We have scanned over the relevant SUSY parameter space, in this case the average 
squark mass rh and gluino mass rrig, assuming SM central values [55]. Moreover, we 
require that the SUSY spectra satisfy the present experimental lower mass bounds [55]. 
In particular, m g > 200 GeV, m > 300 GeV. In addition, we impose that the branching 
ratio (BR) 11 of b — > S7 and the B — B mixing are satisfied at 95% C.L. [59], namely 
2 x 10~ 4 < BR(b -> S7) < 4.5 x 10^ 4 . Then, the allowed ranges for | (5 d LL ) | and | (5f R ) | 

^ ' 23 ^ ' 23 

are obtained by taken into account the above constraints on b — > S7 and B — B mixing. 

In the plots corresponding to QCDF, we have also scanned over the full range of 
the parameters pa,h and <Pa,h m -Y4 and X H , respectively, as defined in Eq.(29). We 

box diagrams in which both up and down mass insertion contribute. See Ref.[ll,51] for further details. 
Then, chargino or gluino contributions cannot be completely disentangled by taking active only one mass 
insertion per time. However, these corrections affect only the Wilson coefficients C\ and C2 in Eq.(7) 
and their effect is quite small in comparison to the other SUSY contributions. In our analysis we can 
safely neglect them. 

11 The branching ratio (BR) of b — > sj is evaluated at the NLO in QCD, as provided in Ref.[50, 56]. 
However, we have not included the 2-loop threshold corrections of SUSY contributions at W scale. For 
these corrections see Ref.[57] for more details. 
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Figure 5: S$k s as a function of arg[(5f jL ) ] (top) and arg[(<5^J ] (bottom) with gluino 

\ / 23 V / 23 

contribution of one mass insertion (8f L ) (top) and (<^r) (down). Left and right plots 

\ / 23 ^ ' 23 

correspond to NF and QCDF, respectively. The region inside the two horizontal lines corresponds 
to the allowed experimental region at 2a level. 
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remind here that Xa and Xh take into account the (unknown) infrared contributions 
in the hard scattering and annihilation diagrams. Regarding the allowed range of p, 
as can be seen from results in Eq.(33), the dominant effect is due to the annihilation 
contributions proportional to Xa- As discussed in section 3, the total width will grow as 
p^log 4 (mj,/A) for large pa • Therefore, as we will show in section 5, by requiring that 
the SUSY contribution to the branching ratio and asymmetries of B — > <pKs is inside the 
la experimental range, an upper bound on pa of order pa — 2 is obtained. 12 

In the corresponding plots evaluated in NF, in order to maximize the SUSY contri- 
butions to CP asymmetry, we have fixed the average of gluon momenta to its minimum 
value (q 2 ) = m 2 / 4. We remind here that {q 2 ) enters as free parameter in matrix element 
of the chromo-magnetic operator Q 8g , see Eq.(23). As we will show later on, the dominant 
SUSY effect to B — > 4>K S amplitude is given by the SUSY contribution to the Wilson 
coefficient C 8g . Therefore, the smaller the (q 2 ) value is, the larger the SUSY contribution 
to the CP asymmetry can be. 

In the framework of NF, the strong phase, which comes from the hadronic matrix 
elements, cannot be predicted. For this reason we set both SM and SUSY strong phases 
to zero. Therefore, the phase 5^ in Eq.(44) can take only the values or 7r, corresponding 
to the relative sign of SM and SUSY amplitudes. 

By comparing the scatter plots in NF and QCDF in Fig. 5 we see that the predictions 
are quite similar. Only the gluino contributions proportional to [5f R j have chances to 



drive SW S toward the region of large and negative values, while the pure effect of (5f 



just approach the negative values region. 

This result can be easily understood by noticing that the dominant SUSY source 
to the B — > <pK s decay amplitude is provided by the chromo-magnetic operator Q 8g . 
In particular, as already mentioned in section 2, gluino contributions to C 8g , which are 
proportional to (J>lr) 23 -> can be very large with respect to the SM ones, being enhanced 
by terms of order mg jm^. In addition, large gluino effects in C 8g may still escape b — > 
57 constraints [58]. This is a remarkable property and it is due to the fact that in the 
gluino contributions to dipole operators Q7 7 ,8 9 , the ratio of \C 8g /Cj 7 \ is enhanced by color 
factors with respect to typical contributions of W or chargino exchanges [48] . 

Now we discuss the chargino effects to S^Kg, which are summarized in Fig. 6. These 
contributions have been first analyzed, in the framework of MIA, in Ref.fll], but only using 
NF approach for evaluating the hadronic matrix elements. Here, we extend our previous 
analysis in [11] by including the corresponding predictions in the QCDF approach. In 
Fig. 6, S^Ks i s plotted versus the argument of the relevant chargino mass insertions, 
namely (8l L ) 32 and (o RL ) 32 . Same conventions as in Fig. 5 have been adopted for left and 
right plots. 

12 We would like to stress here that in the literature the allowed range of pa has been sometime 
overestimated. For instance, in the analysis of [8] the experimental upper bound on BR(B — > <pKs) was 
not imposed, leaving the possibility of larger values of pa ~ 8. 
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Figure 6: As in Fig. 5, but for S$k s as a function of arg[(<5^ L ) 32 ] (top) and arg[(^ L ) 32 ] (bottom) 
with chargino contribution of one mass insertion (5^) 3 2 (top) and (5^) 32 (down). 
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As in the gluino dominated scenario, we have scanned over the relevant SUSY param- 
eter space, in particular, the average squark mass rh, weak gaugino mass M 2 , the \x term, 
and the light right stop mass fhi R . Also tan (5 = 40 has been assumed and we take into ac- 
count the present experimental bounds on SUSY spectra, in particular rh > 300 GeV, the 
lightest chargino mass M x > 90 GeV, and fhi R > 150 GeV. As in the gluino case, we scan 
over the real and imaginary part of the mass insertions (S^ L ) 32 and (Sr L ) 32 , by considering 
the constraints on BR(6 s'j ) and B — B mixing at 95% C.L.. The b — > S7 constraints 
impose stringent bounds on (<5^ L ) 32 , specially at large tan/3 [11]. Finally, as in the other 
plots, we scanned over the QCDF free parameters pa,h < 2 and < 4>a,h < 27r. 

From these results we can see that also for the chargino dominated scenario the pre- 
dictions in NF and in QCDF are quite close, apart from a slight difference in the {5'% L ) Z2 
ones that we will discuss below. The main conclusion in this scenario is that negative 
values of S^Ks cannot be achieved neither in NF nor in QCDF. 

The reason why extensive regions of negative values of S^Ks are excluded here, is 
only due to the b — > constraints [11]. Indeed, as shown in our previous work [11], 
the inclusion of {S1 L ) 32 mass insertion can generate large and negative values of S<pK s , by 
means of chargino contributions to chromo-magnetic operator Q 8g which are enhanced 
by terms of order m x ± /raj,. However, contrary to the gluino scenario, the ratio \Cs g /C^\ 
is not enhanced by color factors and large contributions to Cs g leave unavoidably to the 
breaking of b — > constraints. 

On the other hand, the contribution of (S RL ) 32 is independent of tan/3 and so large 
effects in that could drive S^,k s toward the region of negative values cannot be achieved. 
However, as can be seen from Fig. 6, while both in NF (5^ L ) 32 and {5^ L ) 32 contributions 
are within the 2a experimental range, in NF the (Sr L ) 32 contribution fits better than 
in QCDF. This can be explained by the fact that (S RL ) 32 mainly gives contribution to 
the electroweak operators, whose matrix elements are more sensitive, with respect to the 
other operators, to the approach adopted for their evaluation. 

As shown in our previous work [11], by scanning over the two relevant mass insertions 
(Sr L ) 32 and (<5^ L ) 32 , the b — > constraints on (S^ L ) 32 are a bit more relaxed, but a large 
amount of fine tuning between SUSY parameters is necessary if small values of S<j,K s are 
required. 

In order to understand the behavior of these results, it is useful to look at the nu- 
merical parametrization of the ratio of amplitudes IZ^ = A ^ s IJ(bS<pks) ^ n terms of the 
relevant mass insertions. Below, we present numerical results for lZ<p in both NF and 
QCDF. In QCDF we set to zero the effect of annihilation and hard scattering diagrams, 
corresponding to the choice of pa,h = and 4>a,h = tt- In this case, we expect QCDF 
predictions to be quite close to the NF ones. In particular, for a gluino mass and average 
squark mass of order rh = m g = 500 GeV, we obtain 

~ {-0.08 (5 d LL ) - 120 (S* LR ) } + {L^R}, 
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K,\f DF * {-0.14 x e-^(6 d LL ) 23 - 127 x (5 d LR ) J + {L ^ R} , (55) 

while in the case of chargino, by using gaugino mass M 2 = 200 GeV, \x = 300 GeV, 
mi = 150 GeV, and tan (3 = 40, we have 

n^ F ~ 1.83(^)32- 0-32 (5 U RL ) 32 + 0.41 (5 U LL ) 31 - 0.07 (<^ L ) 31 , 

7^|f^ ~ 1.89 x e- - 07 (^^) 32 - 0.11 x (5 U RL ) 32 

+ 0.43 x e- - 07 (5« L ) 31 - 0.02 x e~^ 17 (5 U RL ) 31 . (56) 

where the first symbol 7Z<f>\^ F means that the corresponding quantity TZ^ has been cal- 
culated in NF approach including only gluino contributions. Analogously for the other 
cases of QCDF and for chargino (x) exchanges. 

From results in Eqs.(55)-(56), it is clear that the largest SUSY effect is provided 
by the gluino and chargino contributions to the chromo-magnetic operator which are 
proportional to (Sf R ^j and (S^ L ) 32 , respectively. However, the b — > S7 constraints play 
a crucial role in this case. For the above SUSY configurations, the b — > 57 decay set the 
following (conservative) constraints on gluino and chargino contributions: | ($lr) 23 1 < 
0.019 and |(5i L ) 32 | < 0.20. Implementing these bounds in Eqs.(55)-(56), we see that 
gluino contribution can easily achieve larger values of = |7^| (see Eqs.(49),(51)) than 
chargino one, and this is the main reason why extensive regions with negative value of 
S^k s are favored and disfavored in Figs. 5 and 6, respectively. 



4.2 CP asymmetry in B — > rj'Ks decay 

Recent measurements of CP asymmetry in B — > r(K s show another discrepancy with SM 
predictions. In particular from Eq.(4), the world average is 

Srf Ks = 0.41 ±0.11, (57) 

and which is about 2.5<r deviation from SM expectations. From these results we see that, 
similarly to what happens in B — > <pKs decay, large deviations from SM are possible. 

Since SUSY contributes to both CP asymmetries with the same CP violating source, it 
is possible that the SUSY effects driving S^Ks towards negative values, could also sizeably 
decrease S v 'k s - The main reason for that is because the leading SUSY contributions to 
the amplitudes of B — > (f>K s and B — > r]'K s enter through the Wilson coefficient of C 8g 
and the operator Q 8g has a comparable matrix elements in both processes. However, since 
NP corrections enter through the quantity ity, the role of the SM contribution will be 
crucial. Indeed, while the B — * <pK s amplitude is purely generated at one-loop in the 
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SM, the B — > rj'Ks one receives tree-level contribution from the SM by means of non- 
vanishing matrix element of Q 2 . Therefore, the increase of SUSY contributions to C$ g is 
now compensated in R v r, by the large SM amplitude contribution. 

We show our results for gluino in Fig. 7, where we have just extended the same analysis 
of B — > 4>K S . Same conventions as in figures for B — > 4>K S have been adopted here. As 
we can see from these results, there is a depletion of the gluino contribution in S v 'k s , 
precisely for the reasons explained above. Regions of negative values of S v 'k s are more 
disfavored with respect to S^Ks, but a minimum of S v >k s — can be easily achieved. 
Comparing NF and QCDF, we can also see that SUSY predictions in NF and QCDF are 
very close. 

Finally, in Fig. 8 we present our results for chargino contributions. Here we see that 
charginos can produce at most a deviation from SM predictions of about ±20 %, and the 
most conspicuous effect is achieved by (5l L ) 32 - These results again show the relevant role 
played by the chromo-magnetic operator. 

As in the case of B — > 4>K S decay, we present below the parametrization of 1Z V > = 
A ^ s S J(B^if Ks) ' k y usm S the same SUSY inputs adopted in Eqs. (56), (55). For gluino 
contributions we have 

*Vir ^ {-0.08 (S d LL ) 23 - 79 (5 d LR ) J - {L^R}, 

n n ,\f DF c {-0.07 x e* a24 (^) 23 - 64(^) 2 J - {L <- R} , (58) 
while for chargino exchanges we obtain 

n^ F ~ 1.2 (5l L ) 32 - 0.01 (S U RL ) 32 + 0.27 (5t L ) 31 - 0.002 (5 U RL ) 31 , 



K V ,\% CDF ~ 0.95 (S U LL ) 32 - 0.025 x e^ ' 19 fe) 32 

+ 0.21 (5 U LL ) 31 - 0.006 x e- i(U9 {8 U RL ) 31 . (59) 

Notice that in the second curly brackets in Eq.(58) there is a minus sign in front. This 
takes into account for the minus sign in the matrix elements of Qi operators as shown in 
Eq.(42). 

In both cases of gluino and chargino contributions, we see that the coefficient of 
(^lr) 23 an d 32 mass insertions is smaller in comparison than the corresponding ones 
in B — > (f)K s , see Eqs. (55), (56). As explained before, this depletion comes from the fact 
that the SM amplitude of A(B — > r/'Ks) receives contribution from tree-level W exchanges 
and it is larger than A(B — > <f)Ks). 

This general behavior is going in the right direction to explain the experimental data. 
In order to show better this effect, we plotted in Figs. 9 the correlations between S^Ks 
versus S v 'k s for both chargino (left plot) and gluino (right plot) in QCDF. For illus- 
trative purposes, in all figures analyzing correlations, we colored the area of the ellipse 
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Figure 7: As in Fig. 5, but for S v >k s - 
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Figure 8: As in Fig. 6, but for S v >k s - 
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Figure 9: Correlation of asymmetries S$k s versus S v >k s with the contribution of one mass 
insertion (<5^) 32 (left) and (Slr)^ (right), for chargino (left) and gluino (right) exchanges. 
Region inside the ellipse corresponds to the allowed experimental ranges at 2a level. 




Figure 10: As in Fig. 9, but for axg[(<^) J=axg[(<^) 23 ] (left) and arg[(<5g H ) 2g ]=arg[(^ £ ) 32 ] 
(right), with the contribution of two mass insertions (<5i#) & (^/J (left) for gluino ex- 

V / 23 V / 23 

changes, and {^lr) & (^11)32 (right) f° r both gluino and chargino exchanges. 
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Figure 11: As in Fig. 9, but taking into account a charged Higgs exchange with mass run = 200 
GeV and tan/3 = 40. 

corresponding to the allowed experimental range at 2a level. 13 In Fig. 9 we scanned over 
the real and imaginary part of the single relevant mass insertions, namely {S^ L ) 32 (left 
plot) and (S^jij (right plot) for chargino and gluino, respectively. 

In conclusion, as can be seen from the results in Fig. 9, pure chargino exchanges have 
no chance to fit data at 2a level, while gluino can fit them quite well. At this point we 
want to stress that sizeable chargino contributions to the CP asymmetries, in particular 
from {Si L ) 32 mass insertion, are ruled out by b — > 57 constraints. Therefore, it would 
be interesting to see if the effects of a light charged Higgs exchange could relax these 
constraints allowing chargino contributions to fit inside the experimental ranges. For 
this purpose, we present in Figs. 11-13 the impact of a light charged Higgs in both chargino 
and gluino exchanges. In particular, plots in Fig. 9 should be directly compared with the 
corresponding ones in Fig. 11, where a charged Higgs with mass mn = 200 GeV and 
tan/5 = 40 has been taken into account. From these results we can see that the effects of 
charged Higgs exchange in the case of 32 mass insertion are negligible, as we expect 
from the fact that terms proportional to (S^ L ) 32 in b — > 57 and b — ► sg amplitudes are not 
enhanced by tan/?. On the other hand, in gluino exchanges with (5lr) 23 or {^ll) 23 ( see 
Figs. 11,12), the most conspicuous effect of charged Higgs contribution is in populating 
the area outside the allowed experimental region. This is clearly due to a destructive 
interference with b — > 57 amplitude, which relaxes the b — > S7 constraints. The most 
relevant effect of a charged Higgs exchange is in the scenario of chargino exchanges with 
(<5^ L ) 32 mass insertion. In this case, as can be seen from Fig. 13, a strong destructive 
interference with b — > 57 amplitude can relax the b — > 57 constraints in the right direction, 

13 All ellipses here have axes of length 4cr. As a first approximation, no correlation between the expec- 
tation values of the two observables have been assumed. 



33 




■1 -0.5 0.5 1 -1 -0.5 0.5 1 

Figure 12: As in Fig. 9, but for gluino contributions with single mass insertion ■ in the 

right plot the effect of a charged Higgs exchange, with mass ran = 200 GeV and tan/3 = 40, 
has been taken into account. 
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Figure 13: As in Fig. 12, but for chargino contribution with single mass insertion ( 5^)32 • 
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allowing chargino predictions to fit inside the experimental region. Moreover, we have 
checked that, for tan/5 = 40, charged Higgs heavier than approximately 600 GeV cannot 
affect the CP asymmetries significantly. 

In a recent paper [9], it has been shown that there exists a particular scenario in which 
gluino contribution can sizeably decrease S^Ks providing on the same time a very modest 
effect in S^Ks- This can be achieved if one assume that both (5lr) 23 and {^rl) 23 mass 
insertions are of the same order, including their CP violating phase. In this case gluino 
contributes with the same weight to both Qi and Qi operators in Eq.(7), where Qi are 
the operators with opposite chirality. Due to the different parity of the final states, as 
explained in section 3, the contribution of the corresponding Wilson coefficients to the 
amplitude of B — > (j)K s and B — > r]'K s enter in combination as Cj + Q and Cj — Q, 
respectively. 

However, the analysis of Ref. [9] was based on the NF approach, and one may won- 
der, if their results still hold in QCDF. For this reason we repeat here the same analysis, 
but in QCDF and scanning over the strong CP phases 4>a,h and pa,h < 2. Results of 
this scenario are shown in Fig. 10 (left plot), where we scanned over two mass insertions 
simultaneously, namely (Sf^j and (5rl) , by assuming that their CP violating phases 
are the same. From these results we can see that a large number of SUSY configuration 
fitting inside the 2a experimental region, can be obtained. 

We have also considered another scenario in which both chargino and gluino exchanges 
are assumed to contribute simultaneously with relevant mass insertions, namely \Sf R ) 



and (Sl L ) 32 . We plot the corresponding results for the correlations between S^k s versus 
SrfKs in the right plot of Fig. 10. As in the case of Fig. 9, we assume a common SUSY CP 
violating phase between the two mass insertions. From these results we can see that also 
in this case a large number of configurations can fit inside the 2a experimental regions. 
This result shows a remarkable fact. The stringent bounds on {S^ L ) 32 from the experi- 
mental limits on BR(B — > X s ^) are relaxed when one considers both gluino and chargino 
contributions, which come with different sign. This does not happen, if two chargino 
mass insertions are considered. Then some configurations with large tan/5 are allowed 
and therefore chargino can contribute significantly to the CP asymmetries S^k s and S v 'k s - 

Finally, we would like to comment about the stability of SUSY predictions for S^k s 
against the low renormalization scale where Wilson coefficients are evaluated. In NF 
the scale dependence in the Wilson coefficients is not compensated by the corresponding 
one in the matrix elements, so a large uncertainty is expected. However, we have noticed 
that also in QCDF this uncertainty still persists. In particular the coefficients in the 
parametrizations in QCDF in Eq. (55)-(56) can vary up to 30-40 % when the fj, b scale is 
changed from /i b — m^/2 to ^ — m h . All the numerical results in this paper correspond to 
the choice ^ = ^1°^ = 4-5 GeV. In both NF and QCDF, this residual scale dependence 
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in is mainly due to the scale dependence in the SM amplitude and in the Cg 9 (/^) 
contribution to the SUSY amplitude. However, we noticed that the main sensitivity to 
the scale in comes from the SM contribution. On the other hand, the main source of 
sensitivity in R n i comes from Csgif^b) in the SUSY amplitude, since the SM one, receiving 
tree-level contributions, is less sensitive to the renormalization scale. 

In the next section we are going to consider the correlations of S^Ks an d SrfKs asym- 
metries versus the corresponding branching ratios of B — > (f)K~s and B — > rj'Ks- 



5 SUSY contributions to BR(B -> <j>(rf)K) 

In this section we discuss the impact of large SUSY contribution, which is required to 
explain the deviations from SM in CP asymmetries of S^m\k s , on the branching ratios 
(BR) of B — > (f)K and B — > rj'K. As shown in Eqs.(5),(6), the experimental measurements 
of these BRs at BaBar, Belle, and CLEO lead to the following averaged results [25]: 



BR(B -> cj)K) = (8.3li;g) x 1(T 6 , 
BR(B -> n'K) = (65.2+^;°) x 1(T 6 . 



(60) 
(61) 



which mean that BR(B — > <f>K) is in good agreement with SM predictions, while the 
experimental value of BR(B — > rj'K) is two to five times larger than the SM one. Since 
these two processes are highly correlated and both are based on the b — > s transition, it 
seems a challenge for the SUSY contribution to enhance BR(B — > rj'K) by a factor of 
two or more, while leaving other BRs and asymmetries inside their experimental ranges. 

The branching ratio of B — > MK, with M = or rj', can be written in terms of the 
corresponding amplitude Mm as 



BR(B -> MK) 



where M M = (MK°\H^ B f =1 \B ) and 



1 \P\ 

8ttMJ 



\M 



1 



M 



Ttot ' 



(62) 



\P\ = 



[(Ml - (m K + m M f) {Ml - (m K - m M ) 2 )] 1/2 



2M B 

The inclusion of SUSY corrections modifies the BR as 



(63) 



BR total (B -> MK) = BR SM (B -> MK) x 1 + 2 cos(^ M - ht)R M + R M 



(64) 



where R 



M 



A SUSY (B^MK) 



A SM (B—> M K) 







M 



arg 



A SUSY (B^MK) 



and 5m is the corresponding strong 



A SM (B->MK) 

phase. The input parameters that we have used in the previous section with pa,h = 1 
and (j)A,H = ir i.e. X a ,h = lead to BR SM (B -> <j>K) = 2.76 x lO" 6 and BR SM (B -> 
V^) = 11.1 x lO" 6 . 
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It is remarkable that in order to produce the experimental result of S^Ks by SUSY 
contribution, the phase 9^ = arg \m(^SjK s f nas ^° °e around ±7r/2 and the strong 
phase 8$ should be of order ir or 0, as shown in section 4. These particular values for the 
phases suppress the leading term in Eq. (64). Therefore, for typical scenarios in which 
R<t, < 1, the SUSY contributions do not enhance BR(B — > 4>K) much. At most a total 
branching ratio of order 3 x BR SM (B — > <f)K) ~ 8.3 x 1CT 6 can be achieved, as we will 
explain below, which is still compatible with the experimental measurements. Clearly, if 
Rff, ^> 1 then the total branching ratio would exceed the experimental bound in Eq.(60). 
However, as shown in section 3, the inclusion of b — > 37 constraints moderates possible 
large SUSY contributions to R^. 

Next we briefly discuss the dependence of the branching ratio on the QCDF free 
parameters Pa,h, related to annihilation and hard scattering diagrams. The only relevant 
parameter here is the pa one, while regarding p H , the BR shows a moderate dependence 
[44]. By using dimensional analysis, one can see that for large values of pa, BR sm (B — > 
(f)K) scales like p\. This result suggests that one can set very strong upper bounds on 
Pa by requiring that BR SM (B — > <f)K) does not exceed experimental range. Clearly, when 
new corrections beyond the SM ones are included, upper bounds on pa could be relaxed 
due to possible negative interferences between SM and new physics amplitudes. 

In order to understand the impact of the annihilation diagrams on BR(B — > (pK), we 
show below the explicit dependence of the SM amplitude on the parameters pa and cf>A'- 

i A SM (B -> <j)Ks) ~ -1.8 x KT 8 + 7.9 x lCT 10 ^ - 6.1 x 10' 10 X 2 A + 2.2 x 10' 10 X H (65) 

where numbers in the right hand side of Eq.(65) are in unity of GeV, and Xa,h are 
defined in Eqs.(29). We have used SM central values as in Table 1 of [44], and omitted 
contributions from imaginary parts which are quite small. 

As can be seen from Eqs.(65) and (29), for pa ^> 1 the SM amplitude could be 
doubled and hence the SM branching ratio would be about four times the result with 
Pa < 1, therefore exceeding the present experimental range. As shown in Ref. [44], and in 
Fig. 14 (see lighter dashed line), the constraint on p A is already obtained for moderate 
values of p A 

It would be interesting to see how much the upper bounds on p A could be modified 
by the inclusion of SUSY corrections, while satisfying all the experimental constraints. 
For this purpose, we plot in Fig. 14 the total branching ratio of B — > Kcj) as a function 
of pa for gluino and chargino contribution. We present the gluino results for the case of 
LR dominant scenario, and the chargino ones with the LL dominant case. The values of 
these mass insertions are varied in the allowed range as done for the other plots. Also 
we scanned over the different values of the strong phases 4>a,h, squark masses, gluino 
and chargino masses. In these plots, the dashed lines correspond to the SM predictions 
for BR(B — > K<p). As can be seen from this figure, by applying the 2a constraints on 
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Figure 14: BR(B — * 4>K) as function of the annihilation parameter pa and for pn = 1, 
<pA,H = 0, for one mass insertion contribution of gluino (sj^ (left) and chargino (5^ L ) 32 
(right). Lighter dashed line corresponds to the SM case. The region inside the two horizontal 
lines corresponds to the allowed experimental range at 2a level. 

BR SM (B — > K(f) one can set very stringent upper bounds on pa, namely p < 2. On the 
other hand, this bound can be relaxed to p < 4 and p < 5, in the case of LL and LR mass 
insertion contributions in up- and down-squark, respectively. However, it is important to 
note that most of the configurations in both gluino and chargino cases that allow for p > 2 
lead to CP asymmetry S^k s outside the 2a range of the experimental results. Therefore, 
as mentioned above, we consider p < 2 as a conservative bound. 

Regarding the BR(B — > rj'K), we find that it is less sensitive to pa than BR(B — > <pK). 
This can be easily observed from the dependence of the amplitude of B — > t]'K on Xa- 
Analogously to Eq.(65), we get 

iA SM {B ->7i'K) ~ 3.6 x 1(T 8 + 0.19 x l(r 10 X A + 4.2x icr 10 x! + 7.7x icr 10 x^ , (66) 

where, as in Eq.(65), numbers in the r.h.s. of Eq.(66) are expressed in units of GeV. Thus, 
comparing Eqs.(65) and (66), we see that the effect of pA in A SM (B — > rj'K) is suppressed 
by two orders of magnitude with respect to the same one in Asm(B — > 4>K). Hence the 
strongest bound one can obtain on pa will come only from BR(B — > <pK). Therefore, in 
all plots of the present work, including the QCDF ones in section 4, we scanned over pa 
by requiring pa < 2. 

As mentioned before, the experimental measurements of BR(B — > rj'K) represent an- 
other large discrepancy with the SM prediction. There have been various efforts to explain 
the large observed branching ratio in the B — > vjK process, based on the peculiarity of 
r)' meson. For instance, intrinsic charm [31] or gluonium contents of rj' [32], have been 
investigated as possible new sources for such an enhancement. Clearly, NP contributions 
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Figure 15: As in Fig.9, but for correlations of S^Ks versus BR(B — > 4>K) (left) and S v 'k s 
versus BR(B — > r/'i^) (right). Gluino contributions with one single mass insertion (^lr) ^ as 
been considered. 

could also be responsible of such discrepancy in BR(B — > i]'K s ), or at least for a part of 
it. In this respect, we are going to analyze next the maximum effect that one can obtain 
from SUSY contributions to BR(B — ► rj'K), by taking into account the experimental 
constraints on the BR(B — ► S^x s , and S v 'k s - 

In Fig. 15 we plot the CP asymmetries S^a^ versus BR(B — > 0if) and versus 
BR(B — > rj'K), where the area of colored ellipse corresponds as usual to the allowed 
experimental region for the correlations at 2a level. We consider here the dominant gluino 
contribution due to (S^r) and scan over the other parameters as before. One can see 
from this figure that, gluino contribution can fit quite well inside the 2a ellipse in both 
left and right plots. It is worth noticing that, gluino effects could also rise BR(B — > <pK) 
to very large values (~ 1.5 — 2 x 10~ 5 ), which are outside the allowed region. 

In the correlation between S v 'k s and BR(B — ► rj'K), with just one mass insertion, 
SUSY contributions can explain the measured large BR(B — > rj'K) and at the same time 
fit the 2a range for S v 'k s - As can be seen from this figure, for S v 'k s within even la 
region, the BR{B — > rj'K) can be large as 7 x 10~ 5 . However, in this case some fine 
tuning between SUSY parameters is necessary. 

In Fig. 16, we present the same correlations as in Fig. 15, but with two simultaneous 
contributions of mass insertions (S^r) and \5rl) in gluino sector. In this case, we can 
easily see that gluino scenario can saturate simultaneously both S^k s and BR(B — > (j>K) 
within their experimental ranges, covering also larger areas with respect to the ones in 
Fig. 15. 

Finally, the combination from gluino and chargino exchanges on BR(B — ► <pK) and 
BR(B — > rj'K) versus the respective CP asymmetries, are shown in Fig. 17. We scan 
on the allowed range of the most relevant mass insertions for these two contributions: 



39 



( §d LFl)23 > ( §d RL)23 



2 4 6 

BR(B->^K s )x10' 



8 10 12 14 
6 



( §d LR)23 > ( §d RL)23 



123456789 10 

^5 



BR(B->ti'K s )x10 3 



Figure 16: As in Fig. 15, but for gluino contributions with two mass insertions {S^j^j a 
ySju^j , with the assumption of arg^tf^^ ] = arg[^ 
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Figure 17: As in Fig. 15, but for gluino and chargino contributions with mass insertions {^'f jR 
and ( s ll) 3 2> with the assumption of axg[[5fj^ ] = arg[(^ L ) 32 ]. 
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Figure 18: Correlations of BR(B — * <pK) versus BR(B — > r/'K) for gluino contributions with 
two mass insertions (<^r) and \S RL ) (left) and for gluino and chargino contributions with 

\ / 23 * 23 

two mass insertions (sf jR j^ and (5^) 32 (right). The 2a constraints on S^Ks and S v ik s are 
applied. 

[Sf R ) (gluino) and (oll) 32 (chargino). We also vary the other parameters as in the 
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previous cases. The message we can learn from these results is that both gluino (LR) 
and chargino (LL) contributions can easily accommodate the experimental results for the 
CP asymmetries and branching ratios of BR(B — > <f)K) and BR(B — > rj'K). As already 
stressed in section 4, this phenomenon is related to the fact that b — > 57 constraints 
are more relaxed for chargino contributions at large tan /3, due to potentially destructive 
interferences of gluino and chargino amplitudes in b — > 37 decay. The main difference 
emerging from results in Figs. 16 and 17, is that for the latter scenario large and negative 
values of asymmetry (S^y ~ —1) can also be easily achieved, even though they are now 
outside the allowed experimental range. 

It is crucial to see, if regions of points which fit inside the two 2a ellipses in Figs. 16 and 
17, actually correspond to the same SUSY configurations. For this purpose, in Fig. 18 
we plot correlations of BR{B — > <pK) versus BR(B — > rj'K) for the same scenarios 
considered in Figs. 16 and 17, where all points satisfy the constraints on S^k s and S^Kg 
at 2a level. As we can see from results in Fig. 18, only the scenario in which both LR 
and RL mass insertions in gluino exchanges contribute, can naturally enhance BR{B — > 
rj'K) inside the allowed experimental range, while respecting all the other constraints 
on CP asymmetries. On the contrary, in the scenario in which both LR(gluino) and 
LL (chargino) are contributing, this enhancement is ruled out by simultaneous constraints 
on CP asymmetries. 

As already mentioned, BR(B — > rj'K) suffers from large theoretical uncertainties due 
to the peculiar structure of rf meson. Therefore, it is not a conservative approach to 
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constrain NP models from the measurements of BR(B — > rj'K), until the role of potential 
new mechanisms responsible of the enhancement of BR(B — > rj'K) in the SM will be not 
completely clarified. 



6 Direct CP asymmetry in b — > 57 versus S^ v ^ Ks 

In this section we analyze the correlation for SUSY predictions between the direct CP 
asymmetry A C p{b — > S7) in 6 — > S7 decay and the other ones in B — > (f)(r}')K s . The CP 
asymmetry in 6 — > 57 is measured in the inclusive radiative decay of 5 — * X s 7 by the 
quantity 

^ - «t) = T(B ^ Xil) + nB ^ Xtl y (er) 

The SM prediction for Acp{b — > 57) is very small, less than 1% in magnitude, but known 
with high precision [46]. Indeed, inclusive decay rates of B mesons are free from large 
theoretical uncertainties since they can be reliably calculated in QCD using the OPE. 
Thus, the observation of sizeable effects in Acp(b 57) would be a clean signal of new 
physics. In particular, large asymmetries are expected in models with enhanced chromo- 
magnetic dipole operators, like for instance supersymmetric models [46]. 

The most recent result reported by BaBar collaboration for A CP (b — > S7) is [60] 

A^T = 0.025 ± 0.050 (stat.) ± 0.015 (syst), (68) 

which corresponds to the following allowed range at 90% confidence level: 

-6% < A b cp sl < 11% . (69) 

Clearly, the present experimental sensitivity is not accurate enough to strongly constrain 
the SM predictions at percent level. 

Recently it has been shown that the SUSY contribution to the CP asymmetry in 
the b — » S7 decay, even with the inclusion of experimental constraints on electric dipole 
moments and branching ratio of B — > X s "f, could be much larger than the SM expectation 
[36] . Therefore, in the light of present experimental results, it would be challenging to 
analyze the SUSY predictions for the correlation between A C p{b — > 57) and S^^Ks, 
since in SUSY models these asymmetries are strongly correlated. 

The relevant operators of the effective Hamiltonian in Eq.(7) that play a role in 
Acp{b — > 57), are given by Q2, Qi^iQsg- We remind here that Q2, defined in Eqs.(8), 
is the usual current-current operator induced at tree-level. Then, the expression for 
Acpip — > S7) at the NLO accuracy, is given by [46] 
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|CV 7 | 2 [81" 



A CP (b^ S1 ) = -^f\—lmC 2 Cl- — [v(z) + b(z,5)]lm(l + e s )C 2 C^ 
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-b(z,5)Im[(l + e s )C 2 C* g \ , (70) 
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where z = (m c /m b ) 2 and v(z) is [46] 



v(z) = (^ + \nz + \n 2 z-^Pj + (^\n 2 z - z + - Unz^j z 2 + 0(z 3 ). (71) 

The function b(z, 5) is defined as b(z, 5) = g(z, 1) — g(z, 1 — 5), where the parameter 5 
is related to the experimental cut on the photon energy, > (1 — 5)m&/2, and finally 
g(z,y) is given by 

3y(y - 2z 



g(z, y) = 6(y - Az) j (y 2 - Ayz + Qz 2 ) In (JI- + ^-1 




In the SM, the Wilson coefficients C 2 , C 77 , and C 8s are real, in particular Qf M ~ 1.08, 
C 77 = —0.318, and C Sg = —0.151, so that the only source of CP violation comes from the 
e s parameter inside Eq.(70) which is defined in terms of the CKM matrix elements as 

= W^A 2 (^-p)~O(10- 2 ). (73) 

Vts V tb 

In general SUSY models, CV 7 and C$ g may be complex, and the corresponding phase 
can provide the dominant CP violating source to A 1 ^ 1 . In MIA framework, the SUSY 
contributions to the Wilson coefficients C 7l and C Sg are given in terms of (^5)23 (for 
gluino) and (<5^ B ) 32 (for chargino). Therefore, the SUSY phase in these mass insertions 
is the only source of CP violation in this process. As pointed out in previous sections, in 
general SUSY models, effects induced by the dipole operators Q77 and Q$ g , which have 
opposite chirality to Q77 and Q$ g , are non-negligible. In the MSSM these contributions are 
suppressed by terms of order 0{m s jm^) due to the universality of the A-terms. However, 
in general models one should take them into account, since they might be of the same 
order than Cy 7 and C 8g . Denoting by C 7j and C Sg the Wilson coefficients multiplying the 
new operators Q77 and Qs g , the expression for the asymmetry in Eq.(70) will be modified 
by making the replacement 

dC* ^C t C* + C t C*. (74) 

Notice that, since 14 C 2 = 0, the only modification in the numerator of Eq.(70) is 
due to the new term C-j^C^j- However, if only one single mass insertion is taken into 
account, both CV 7 and C$ g will be proportional to the same mass insertion, and therefore 
Im (CggCj^j = 0. Therefore, the effects of these new operators will enter only through 
C7 7 in the denominator, by means of the shift |C7 7 | 2 — > |C7 7 | 2 + |C7 7 | 2 . 

It is worth mentioning that in the scenarios dominated by one single mass insertion, 
the leading SUSY contribution to the CP asymmetry A C p(b — > S7) is due to the C 77 . This 

14 In principle, C2 might receive some contributions from RGE due to the mixing of the other Qi 
operators with Q2. However, the radiative effects induced by C\ on C2 are quite small and so we will 
neglect them in our analysis. 
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can be understood as follows. Using our previous inputs for charm and bottom masses and 
also assuming, for instance, an energy resolution 5 ~ 0.3 so that b(z,S) ~ 0.17, Eq.(70) 
can be written as 

A CP {b - s 7 ) ~ (0.045 Im[C 7 * 7 ] - 0.44 lm[C 7l C* g ] + 0.006 Im[C* g ]) . (75) 

We remind here that the low energy coefficients CV 7 and C£ g at //& scale are related to the 
high energy ones by 

C 7 S , USY (m b ) * 0.66 C™"W) + 0.09 (76) 
Cl^im,) * 0.69 C7™*W). (77) 

As shown in section 2, the kind of mass insertions that will appear in CV 7 (/iw) are the 
same that determine Cg g (nw)- I n the case of SUSY contribution from one single mass 
insertion, both CV 7 (/^) and (7| fl (//&) will acquire the same phase and so the second term 
in Eq.(70) identically vanishes. 

In conclusion, the relevant mass insertions contributing to A C p(b — > s 7 ) give also 
the dominant effects in S^Ks an d S v 'k s - Therefore, it is expected that these three CP 
asymmetries should be highly correlated, since all of them depend on the same source of 
SUSY CP violation. 

In Fig. 19 we plot the correlations between Sk s 4> and A C p{b — > 57) and also S v 'k s 
and Acp(b — > 57) for gluino exchanges with one mass insertion, namely (5lr) 23 - As 
we can see from results in Fig. 19, a specific trend emerges from this scenario. The 
simultaneous combination of S^Ks and S v >k s constraints, favors the SUSY predictions to 
be in the region of positive values of Acp(b — > S7) asymmetry. 15 As shown in Fig. 20, 
analogous results are also obtained when two mass insertions in gluino sector are taken 
simultaneously, namely (5lr) 23 an d ($lr) 23 - However, in this case also negative values 
of A b cp~' are allowed. Although these points are favored by S,p Ks , they are disfavored by 

SrfKs- 

Finally, in Fig. 21 we show our results for two mass insertions (5f^j and (5^ L ) 32 
with both gluino and chargino exchanges. In this case we see that S^Ks constraints do 
not set any restriction on Acp(b — > s 7 ), and also large and positive values of Acp{b — > s 7 ) 
asymmetry can be achieved. However, by imposing the constraints on S v /k s , see plot on 
the right side of Fig. 21, the region of negative Acp(b — > s 7 ) is disfavored in this scenario 
as well. 



15 This result is also in agreement with the corresponding predictions for the correlation of S^Ks versus 
Acp (b — > S7) in the article of Kane et al. in Ref. [6] . 
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Figure 19: Correlations of S^Ks versus Acp(b — > 57) (left) and S^Ks versus j4cp(6 — ► S7) 
(right), for gluino contributions with one single mass insertion ySf j j^j i 
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Figure 20: As in Fig. 19, but for gluino contributions with two mass insertions (^lr) 23 anc ^ 
(^rl) 23 - 
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Figure 21: As in Fig. 19, but for gluino and chargino contributions with mass insertions (<$£#j 
and (Sll) 32 respectively. 

7 Conclusions 

The 5-meson decays to (f>K, rj'K, and to X57 provide a clean window to the physics 
beyond the SM. In this paper we have systematically analyzed the supersymmetric con- 
tributions to the CP asymmetries and the branching ratios of these processes in a model 
independent way. The relevant SUSY contributions in the b — > s transitions, namely 
chargino and gluino exchanges in box and penguin diagrams, have been considered by 
using the mass insertion method. 

We have provided analytical expressions for all the relevant Wilson coefficients. Both 
naive factorization and QCD factorization approximation to determine the hadronic ma- 
trix elements have been employed. We showed that the SUSY predictions for the mixing 
CP asymmetry of S^Ks an d S v >k s are not very sensitive to the kind of approximation 
adopted for evaluating hadronic matrix elements. 

We found that due to the stringent constraints from the experimental measurements 
of BR(b — > 57) , the scenario with pure chargino exchanges cannot give large and negative 
values for CP asymmetry S^,k s - Indeed, by combining present experimental constraints 
on S^Ks an d SrfKs asymmetries at 2a level, pure chargino contributions can be almost 
ruled out. On the other hand, it is quite possible for gluino exchanges to account for 
S^Ks and SrfKs a t tli e same time. Interestingly, we have shown that the charged Higgs, 
if not too heavy, may change the chargino and gluino contributions to enhance the CP 
asymmetries considerably. 

The branching ratios of these decays have also been considered. It has been noticed 
that the p parameter is strictly bounded by the B — > <pK branching ratio, and this 
influences strongly the possible SUSY contributions to the asymmetries. Investigating 
the correlations between CP asymmetries and branching ratios, pointed out that the 
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most favored SUSY scenarios, which easily satisfy all the experimental results, are the 
ones with inclusion of two mass insertions. In particular, a pure gluino dominated scenario, 
in which both LR and RL mass insertions are contributing and the one in which both 
gluino and chargino contribute with LR and LL mass insertions in down- and up-squark 
sectors, respectively. In the latter scenario we show that chargino exchanges provide 
sizeable effects. This is due to the fact that b — > S7 constraints could be relaxed by 
potentially destructive interferences between gluino and chargino amplitudes. Finally, it 
is remarkable to notice that in the scenario in which both LR and RL down squark mass 
insertions dominate, the observed large enhancement of BR(B — > rj'K) could be naturally 
explained, while respecting all the other experimental constraints on CP asymmetries. 

We also discussed the correlations between the CP asymmetries of these processes and 
the direct CP asymmetry in b — > S7 decay. In this case, we found that the general trend 
of SUSY models, satisfying all the experimental constraints, strongly favors large and 
positive contributions to b — > 57 asymmetry. 

More precise measurements would certainly allow us to draw more definite conclusions 
on the SUSY models that can accommodate these data. 



Acknowledgments 

The authors thank Debrupa Chakraverty for her efforts in the initial stages of this work. 
We would like to thank G. D'Agostini, E. Kou, and K. Osterberg for very useful dis- 
cussions. E.G. acknowledges the CERN PH-TH division and IPPP, CPT of Durham 
University for their kind hospitality during the preparation of this work. The authors 
appreciate the financial support from the Academy of Finland (project numbers 104368 
and 54023) and from the Associate Scheme of ICTP. 



47 



Appendix 

A Chargino contributions in MIA 

Here we provide the expressions for chargino contributions, at leading order in MIA, to 
the R F B quantities in Eq.(15), where F refers to D,E,C,B^ d ), and M 7)9 [11]. 
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R B L d = E vzyiMtfxwjPiixjtXv) 



i,J=l,2 
i,j=l,2 

= El^il'^^.^-nE^i^^^p^^) 

777- 

i i m b 

R M^, g = ~ E V tt V i2 X Wi P M^, g ( X i) 

i 

R m*, 9 = EW*™^.><)- (78) 

i 

Regarding the R F terms at the zero order in mass insertion, we have 
R° D = 2x w £ \V t2 \ 2 D x (x t ) 

i=l,2 



i=l,2 

r° c = 2 x: c/ac/;iw; 2 cf(^,^) 

Z ij=l,2 

= E x Wj V i2 V* 2 Ut l U* 1 B^\x j ,x j ,x ij ) 
r m, = -xw E \V l2 \ 2 [fM) + \F2(x^ 

r m 9 = -xw E \y^\ 2 F2{xi) (79) 

i=l,2 

where Y b is the Yukawa coupling of bottom quark, x wi = m^/m 2 ,^ x w = m^/m 2 , 
Xi = m^Jrh 2 , Xi = m 2 /m 2 ., and x^- = ro^./m^.. The chargino mixing matrices {7 
and V, are defined in Eq.(13). The loop functions of penguin Pd,e,c, ^xi^xi^xi b° x 
B^ u ' d ^\ pjv> d >9) ^ anc j ma gnetic and chromo-magnetic penguin diagrams iq,2, Pj^ , -Fm^ 
are reported in appendix C. 16 



16 The minus sign appearing in front of the right-hand-side of equation above for R% d , takes into account 
the correction of a sign mistake in the first reference of [51] for the chargino contributions to down- type 
box diagrams. This mistake was also pointed out in the second reference of [51]. This sign correction has 
also been re- absorbed in the P% function given in appendix C. 
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B Chargino contributions in MIA for the light Ir case 

In this appendix we provide the relevant formulas for the case in which the mass of the 
stop-right ) is lighter than the average squark mass fn 2 . We recall here that this will 
modify only the expressions of R RL , R RR , and R%, since the light stop right does not affect 
Rp L . In the case of R RR E c and R° D E c the functional forms of R RR E c and R° D E c remain 
unchanged, while the arguments of the functions involved are changed as xyy — > xwt, 
Xi Xit, and Xi — > x^, where xwt = m w/ m f R ^ x u = m 'x l l m \ R i xu = m i R / m xi' Regarding 
the box-type contributions 17 RB R d an d R% ud -> the functions P^ ,d ^ and B x u,d ^ appearing 
there should be replaced as follows 



,(it,d) 



(Xj , Xij ) 



{Xjt, Xj, Xij) 



B^ ' ^ (Xj, Xj, Xij) > B^ ' ^ (x jt, X j , Xij) , 



(80) 



where 



P^(x,y,z) 



l : x^-B^ d \x,y,z) 



2~dx~* — (81) 
In the case of R F R and R RL the analytical expression of loop functions of penguin Pd,e,c, 
box P B u,d \ and magnetic and chromo-magnetic penguin diagrams Pj^ g and P^ R respec- 
tively, should be replaced with the following ones 



Pd,e( x i) 



P$' 2 \xi,Xj) 



PD,E{Xi,Xit) 

Pci^i, Xit) 

p(l>2)/ \ 



(82) 



where 



P D (Xi,Xit) 

Pe{x%i Xu) 



*C \Xii Xu, Xj , Xjt) 



Pn ^ (Xi, Xit) 



C 

>(«)f~ ~ 



TD\^) ( ry> ry> ry> \ 

1 B K^J i Jj 3ti ""13 ) 



_ y [xi t D x (x it ) - XiD x (xi)} 
2 

[x it E x (xit) - XiE x (xi)] 



(x t - 1) 
2 

S - 1) 
4 

fa " 1) 
1 

2(x t - 1) 
1 



2(x t - 1 



' ^ {Xjti Xu) ' H^'j'^'i) 

) -^'(Si,^)] 

Xj, Xij) — B x \xj, Xj, Xij) 
B x ^ {.Xjt, Xj, x^) — B x •* (xj, Xj, Xij 



17 Here we have corrected the formula in [11]. The effects are 0(X 2 ) and do not affect the numerical 
results. 
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^ > M 1 ( X ii X it) 
^(^i) x it) 



(x t - 


1) 


1 




(x t - 


1) 


1 




(x t - 


1) 


1 




(x t - 


1) 



x it (f^xu) + -F 2 (x it )j - x { (Fifa) + -F 2 (xi) 
xu (F 3 (x lt ) + 2 -F i (x lt ))-x l (F s (xO + f F 4 (x 4 ) 



Pjfi(xi,x a ) = J2 - ^ [x it F 2 (x it ) - x i F 2 (x i )\ 

P^{xi,x it ) = ^ [x it F A (x it ) - XiF 4 (xi)] . (83) 

The expressions for the loop functions involved in Eqs.(81) and (83) can be found in 
Appendix C. 

C Basic SM and SUSY Loop functions 

Here we report the loop functions entering in the SM and SUSY diagrams at 1-loop. 
• SM functions 



x 2 25 - 19x) -3a; 4 + 30a; 3 - 54a; 2 + 32a; - 8 , 

B(x) = — ^ ^ + - J -log a; 

36 (a;-!) 3 18 (x - if 



x (x - 6) x (3x + 2) 

= ,(^ + 11,-18) + 
k ' 12(i-l) 3 6(i'-l) 4 

BW = -W^T) + W^Tf-' 0ix (84) 



Charged Higgs functions 



^ , N x 47a; 2 - 79x + 38 x -3a; 2 + 6x - 4 ) , 

D H U = — -o -H -j — -log a; 

108 (x -l) 3 18 (x -l) 4 

C H (x) = -^B(x) 

^ , , a; (7x 2 - 29a; + 16) x (3x - 2) n 
36 (a; — 1) 6 (x — 1) 
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• Chargino functions in MIA 



P D (x) 
Pe(x) 

Pb(x,v) 



+ 



+ 



Pm*{*) 



2x (-22 + 60x- 45x 2 + 4x 3 + 3x 4 -3 (3-9x 2 + 4x 3 ) logx) 

27 (1-x) 5 

x (-l + 6x-18x 2 + 10x 3 + 3x 4 -12x 3 log a;) 

9(1 -xf 
x (3 — 4 x + x 2 + 2 log x) 



1 



8 (x - y) 
4 (x - y) 



(1-x) 3 

x 2 (x - 1 - logx) _ y 2 (y - 1 - logy) 

(x - l) 2 (y - l) 2 

x (x — 1 — log x) y (y — 1 — log y) 



(x-l) ; 



(l/-l) J 



-y — x (1 — 3 x + y) x (x 3 + y — 3 xy + y 2 ) log a; 



4(x-l) 2 (x-y) 2 



2(x-l) 3 (x-y) 3 



xy logy 



2(x-y) 3 (y-1) 
x (3 y - x (1 + x + y)) x (x 3 + (x - 3) x 2 y + y 2 ) logx 



4(x-l) 2 (x-y) 2 
xy 2 logy 



2(x-l) 3 (x-y) 3 



2(x-y) 3 (y-1) 

x (-2 - 9x + 18x 2 - 7x 3 + 3x (x 2 - 3) logx) 
9(1 -x) 5 

x (13 - 20x + 7x 2 + (6 + 4x - 4x 2 ) logx) 
6 (1-x) 4 

x (-1 + 9x + 9x 2 - 17x 3 + 6x 2 (3 + x) logx) 

12(1 -x) 5 
x (—1 — 4x + 5x 2 — 2x (2 + x) logx) 
2 (1-x) 4 



(86) 



• Chargino functions in MIA with a light Ir 



D x (x) 



-43x 2 + lOlx - 52) (2x 3 - 9x + 6) , 



108 (x - 1) 



16 (y - x) 



18 (x- 1) 

x 2 y 2 

log x log y 



x- 1 
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E x (x 



B^(x,y,z 
B^(x,y,z 

F(x,y,z 
F 1 (x 
F 2 {x 
F 3 (x 
F 4 (x 



xy 



S(y-x) 



x 



x — 1 



log a; 



-11a; 2 + 7a; - 2) 



+ 



X 



logy 



log a; 



36 (x - l) 3 6 (x - 1) 

- F(x, y, z) + F(y, z, x) + F(z, x, y) 

xF(x, y, z) + yF(y, z, x) + zF(z, x, y) 

x log a; 
(x — l)(x — y)(x — z) 
(x 3 - Qx 2 + 3x + 2 + Qx \ogx) 

12 {x- l) 4 
(2a; 3 + 3a; 2 - 6x + 1 - 6a; 2 log a;) 

12 (a; -l) 4 
(a; 2 -4x + 3 + 21oga;) 

2(x- l) 3 
(a; 2 — 1 + 2a; log a;) 
2(a;-l) 3 



Gluino functions in MIA 



Pi(x) 
P 2 (x) 
Mi(x) 
M 2 (x) 
M 3 (x) 
M 4 (x) 



1 - 6x+ 18x 2 - 10a; 3 - 3x 4 + 12a; 3 ln(a;) 
18 (x - l) 5 

7- 18a; + 9a; 2 + 2a; 3 + 3 ln(x) - 9x 2 ln(x) 

9(x-lf 
1 + Ax - 5x 2 + Ax \n(x) + 2x 2 \n(x) 

2(1 - a;) 4 
-5 + Ax + x 2 - 2 \r(x) - Ax \n(x) 
2(1 -x) 4 

-1 + 9a; + 9a; 2 - 17a; 3 + 18a; 2 ln(x) + 6a; 3 \n(x) 
12 (x- l) 5 

— 1 — 9 a; + 9 a; 2 + a; 3 — 6a; ln(a;) — 6a; 2 ln(a;) 



6 (x - if 



Bi(x) = ^Mi(i) 
B 2 (x) = -xM 2 (x) 
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